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JU  SUMMARY 

An  experimental  investigation  has  been  carried  out  to  determine  the 
effect  of  heating  on  boundary  layer  transition  in  a liquid.  For  this  purpose 
water  was  circulated  through  a steam-heated  smooth  brass  tube  over  a range  of 
diameter  Reynolds  numbers  from  30,000  to  140,000.  With  the  disturbance  level 
present  in  the  system  it  was  possible  to  produce  adiabatic  length  Reynolds 
numbers  of  transition  up  to  850,000.  The  location  of  transition  was 
determined  from  measurements  of  local  apparent  friction  factor  along  the 
tube,  and  also  by  the  use  of  a boundarv  layer  probe  which  detected  the 
change  in  velocity  profile  at  the  onset  of  transition.  Heating  was  found  to 
have  no  significant  effect  on  the  position  of  transition,  and  resulted  in  a 
slight  decrease  of  free-stream  length  Reynolds  number  of  transition*  while 
the  mean  boundary  layer  Reynolds  number  of  transition  was  slightly  increased. 
Small  roughnesses  on  the  boundary  walls  were  found  to  have  a very  pronounced 
effect  on  the  location  of  transition. 

The  results  indicate  that  for  an  ordinary  water  flow  system,  the  use  of 
heating  could  not  be  expected  to  appreciably  cause  stabilization  of  the 
laminar  boundary  layer. 

An  expression  for  local  apparent  friction  factor  in  the  laminar  entry 
of  a tube  has  been  developed  by  the  use  of  an  approximate  boundary  layer 
computation  method.  The  relation  derived  applies  within  a range  where 
Re^Rej^  is  less  than  10“  . 

A second  solution  for  apparent  friction  factor,  which  applies  throughout 
the  whole  laminar  entry  region  of  a tube,  has  been  derived  by  substituting  an 
approximate  velocity  profile  shape  into  the  continuity  and  momentum  equations 
The  velocity  profile  chosen  approached  flat  plate  behavior  near  thejtube 
entrance,  and  gradually  changed  toward  a parabolic  shape  as  the  flow 
developed.  Graphical  integration  was  used  to  obtain  the  final  results  in 
tabular  form. 
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II.  INTRODUCTION 

Brief  Introduction  to  Mathematical  Stab 11 ltv  Theory 

The  current  concept  of  the  process  governing  the  transition  of  a 
boundary  laver  from  laminar  to  turbulent  motion  is  based  on  stability 
considerations . At  a certain  critical  value  of  the  Reynolds  number,  the 
flow  becomes  unstable  for  small  disturbances,  which  then  amplify  until  the 
motion  becomes  turbulent.  The  point  at  which  the  critical  Revnolds  number 
is  attained  is  termed  the  instabilitv  point,  while  the  point  at  which  the 
boundary  layer  begins  to  deviate  from  laminar  motion  is  called  the  transition 
point.  Between  the  two  lies  the  zone  in  which  small  disturbances  undergo 
amplification.  Since  the  mathematical  investigations  deal  with  the  criteria 
for  the  growth  of  a small  disturbance  in  the  stream,  it  should  be  noted 
that  they  predict  the  location  of  the  instability  point,  and  not  that  of 
the  transition  point. 

The  general  approach  to  the  problem  has  been  to  superimpose  a small 
disturbance  motion  upon  the  differential  equations  governing  the  basic  flow. 
Only  first  order  terms  of  the  small  perturbation  velocity  are  retained,  and 
the  amplitude  of  the  disturbance  motion  is  investigated  as  a function  of 
time.  These  calculations,  carried  out  by  Tollmien^ 23)«<-,  schlichting^7^ , 
and  others,  have  shown  that  the  critical  Reynolds  number  is  a function  of 
the  frequency  of  the  disturbance,  and  also  of  the  shape  of  the  velocity 
profile.  If  the  behavior  for  profiles  with  and  without  inflection  points 
is  compared,  it  is  found,  in  general,  that  those  having  inflections  are 
unstable  over  a wider  range  of  frequencies  than  those  without  inflections. 
Moreover,  the  value  of  the  critical  Reynolds  number,  below  which  the  flow 
is  completely  stable  for  all  disturbances,  is  lower  in  the  former  case. 

Since  the  shape  of  the  velocity  profile  can  be  altered  by  pressure  gradients 
in  the  direction  of  motion,  or  by  temperature  distribution  across  the 

# Superscripts  in  parenthesis  refer  to  reference  numbers  in  the 
bibliography. 
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boundary  layer,  these  factors  have  an  effect  on  the  predicted  stability  of 
the  flow , 

Usinrr  the  extension  of  the  theory  for  the  case  of  a compressible 
fluid^*^  , numerical  calculations  have  been  carried  out  for  gases  flowing 
past  a flat  plate.  The  calculations  predict,  in  the  case  of  air,  an 
increase  in  stability  when  the  plate  is  cooled  and  a decrease  when  the 
plate  is  heated.  This  comes  about  because  an  inflection  point  tends  to  be 
developed  in  the  velocitv  profile  for  the  case  of  a gas  flowing  over  a hot 
plate.  Since  the  rate  of  change  of  viscosity  with  temperature  in  a liquid 
is  negative,  whereas  in  a gas  it  is  positive,  it  micht  be  expected  by 
analogy  that  the  heating  of  a plate  with  liquid  flowing  over  it  would  tend 
to  increase  laminar  stability,  although  this  case  has  not  been  adequately 
treated  at  present. 

Summary  of  Experimental  Work 

(19) 

Adiabatic  experiments  carried  out  by  Schubauer  and  Skramstad'  ' have 
fully  substantiated  the  theoretical  dependence  of  stability  on  both  the 
frequency  of  the  disturbance,  and  the  external  pressure  gradient.  Experi- 
ments by  Liepmann  and  Fila^10^  have  shown  that  the  boundary  layer  velocity 
profile  for  air  passing  over  a heated  flat  plate  develops  an  inflection, 
and  that  the  transition  point  is  caused  to  move  upstream  as  heat  is  applied. 
Experiments  by  Scherrer^15^  for  cooling  an  air  stream  flowing  externally 
over  a cone  at  Mach  number  greater  than  unity,  have  shown  a small  Increase 
in  the  transition  Reynolds  number,  but  the  increases  were  much  less  than 
that  predicted  by  theory.  In  the  case  of  internal  flows,  recent 
investigations  by  Kline  and  Shapiro^6)  for  cooling  air  flowing  within  a 
smooth  tube  have  shown  no  appreciable  effect  on  the  transition  point  even 
when  th^  wall  temperature  was  reduced  over  200°F  below  the  stream  air 
temperature.  These  results  clearly  demonstrate  that  more  experimental 
work  on  the  effect  of  temperature  gradients  on  stability  is  needed  before 

a full  understanding  can  be  obtained. 
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Nature  of  the  Present  Investigation 
In  the  present  work,  water  was  circulated  through  a tube  with  heated 
walls  to  determine  whether  transition  could  be  delayed  by  this  means.  If 
successful,  this  technique  would  have  application  in  reducing  the  surface 
drag  on  under-sea  projectiles,  on  submarines,  and  in  the  inlets  of 
hydraulic  machines.  Since  the  mathematical  theory  predicts  only  the 
instability  point  for  motion  with  vanishingly  small  disturbances,  it  was 
necessarv  to  resort  to  experiment  to  determine  whether  the  theoretical 
trends  for  stability  apply  also  to  the  behavior  of  the  transition  point 
in  a fluid  stream  with  a small  but  finite  disturbance  level. 
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III.  DESCRIPTION  OF  THE  TEST  APPARATUS 
Flow  System 

A schematic  representation  of  the  water  recirculating  system  used  for 
this  investigation  is  shown  in  Figure  1.  Water  drawn  from  the  city  main 
was  contained  within  a five  hundred  gallon  storage  tank  and  circulated 
through  the  test  equipment  by  means  of  a five  horsepower  centrifugal  pump. 

A by-pass  back  into  the  tank  was  provided  at  the  pump  discharge  for  the 
purpose  of  flow  regulation.  The  piping  was  arranged  so  that  the  entire 
flow  could  be  diverted  at  low  velocity  through  a filter  to  be  carefully 
cleaned  before  introduction  into  the  test  section.  A standard  orifice  meter 
was  located  in  the  main  line  and  mounted  according  to  ASME  specifications'  ' 
Stilling  Chamber 

From  the  main  line  the  flow  passed  through  a conical  diffuser  fabricated 
of  galvanized  sheet  metal.  The  diffuser  helped  to  prevent  the  formation  of 
large  scale  turbulence  during  the  transition  from  the  three  inch  pipe  line 
to  the  ten  inch  diameter  stilling  chamber.  Following  the  diffuser  the  water 
entered  a ten  inch  pipe  tightly  packed  with  approximately  two  thousand  thin 
walled,  extruded  aluminum  tubes  having  a length  diameter  ratio  of  forty 
(Figure  2).  This  honeycomb  was  sized  to  insure  laminar  flow  in  and  between 
the  tubes  for  all  flow  rates  that  were  anticipated  during  the  tests,  and 
was  employed  to  remove  any  swirling  motion  from  the  fluid. 

A stilling  chamber  containing  a series  of  fifteen  turbulence  damping 
screens  was  located  after  the  honeycomb.  The  chamber  design  was  such  that 
it  should  theoretically  provide  about  a two  hundred  fold  reduction  of  the 
entering  turbulence  intensity.  The  screens  chosen,  according  to  the  data 
of  references  2,  3,  4,  and  20,  were  34  mesh  and  had  a wire  diameter  of 
0.0065  inches.  Each  was  tightly  stretched  and  soldered  to  an  interlocking 
annular  ring  as  shown  in  Figures  2 and  3.  The  screens  were  spaced  two 
inches  apart  so  that  any  low  scale  turbulence  generated  by  a screen  would 
be  damped  before  the  next  screen  was  reached.  The  combination  of  chamber 
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size,  wire  diameter,  and  mesh  number,  was  chosen  such  that  a critical 
Reynolds  number  would  not  be  exceeded,  that  is,  vortices  would  not  be  shed 
from  the  screen  wires  for  the  range  of  flow  rates  used.  Before  entrance 
into  the  test  section,  a further  reduction  in  turbulence  was  obtained  due 
to  the  large  reduction  in  cross-sectional  area  when  the  stream  was  contracted 
from  the  ten  inch  diameter  chamber  to  the  two  inch  diameter  test  section. 
Contraction  Nozzles 

Following  the  stilling  chamber  the  flow  passed  through  a double 
contraction  nozzle,  going  first  from  10  to  4 1/2  inches  and  then  from 
4 1/2  to  2 inches  in  diameter.  Each  contraction  was  carefully  machined  to 
a double  cubic  contour  according  to  the  data  of  Rouse  and  Hassan^4^,  so  as 
to  avoid  any  unfavorable  pressure  gradients  which  might  induce  boundary 
layer  separation.  The  finish  machining  was  completed  with  the  two  contrac- 
tions clamped  together  so  that  a smooth  joint  would  be  assured,  and  the 
whole  interior  was  buffed  to  a high  luster. 

The  wall  of  the  second  contraction  was  hollow  (Figure  3),  and  cooling 
water  co\ild  be  circulated  through  it  during  the  heated  runs  in  order  that  it 
might  be  maintained  at  the  same  temperature  as  the  water  entering  the  test 
section.  This  was  necessary  because  the  contraction  could  not  be  effectively 
thermally  isolated  from  the  steam  jacket  used  to  heat  the  test  pipe.  With- 
out cooling  of  the  contraction  walls,  there  would  have  resulted  an  undesir- 
able thermal  boundary  layer  formation  in  the  stream  before  reaching  the 
test  section. 

Design  of  the  Test  Section 

The  test  section  was  constructed  from  a 2 inch  internal  diameter  drawn, 
seamless  brass  tube  96  inches  long  and  with  a 0.125  inch  wall  thickness. 

At  the  entrance  of  the  section  a boundary  layer  suction  slot  was  provided, 
through  which  up  to  10  percent  of  the  flow  could  be  removed.  The  slot  was 
designed  according  to  specifications  obtained  from  reference  12,  and  the 
configuration  used  is  shown  in  Figure  4.  The  water  removed  passed  from  the 
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annular  collection  ring  surrounding  the  slot,  through  eight  copper  tubes  of 
equal  length,  and  thence  through  the  end-plate  of  the  steam  jacket  into  a 
collection  header.  The  boundary  layer  flow  was  metered  by  a standard 
flange -tap  orifice  plate,  and  was  returned  to  the  main  storage  tank  through 
an  appropriate  piping  system  and  control  valve. 

The  test  pipe  was  provided  with  18  wall  static  pressure  taps,  each 
having  a diameter  of  0.020  inches  (Figure  4).  These  were  spaced  at 
distances  in  proportion  to  their  distance  from  the  boundary  suction  lip, 
except  for  the  first  few  taps  which  were  uniformly  spaced  2 inches  apart 
(see  Table  I).  To  prevent  the  build  up  of  a disturbance  from  one  tap  to 
the  next,  the  taps  were  successively  spaced  from  one  another  at  a helix 
angle  of  110  degrees,  so  that  no  two  taps  were  in  line  with  the  axial 
direction.  The  taps  were  carefully  machined  by  alternate  drilling  and 
lapping  operations  to  make  them  as  uniform  and  square  as  possible,  and  to 
eliminate  burrs  at  the  inside  edges.  After  the  final  drilling  operation  the 
whole  tube  was  lapped  and  given  a high  polish  on  the  inside. 

Six  thermocouples  were  mounted  at  intervals  in  tangential  slots  1/16 
inch  deep,  spaced  on  alternate  sides  of  the  test  pipe  (see  Table  II).  The 
thermocouples  were  placed  in  the  slots,  and  the  slots  filled  with  solder 
which  was  filed  down  flush  with  the  outer  tube  wall. 

The  whole  test  section  was  enclosed  within  a steam  jacket  made  of 
standard  eight  inch  pipe.  The  jacket  was  provided  with  steam  at  3 psi  gage. 
All  of  the  pressure  tap  and  thermocouple  connections  were  brought  out 
through  the  steam  jacket  end-plate.  The  exit  end  of  the  test  section  was 
fixed  to  the  steam  jacket  through  a movable  0-ring  seal.  This  was  done  to 
prevent  thermal  stresses  arising  from  the  unequal  expansion  of  the  jacket 
and  test  pipe. 

All  of  the  return  outlets  in  the  storage  tank  were  extended  below  the 
water  surface  to  prevent  splashing  and  subsequent  entrainment  of  air  bubbles 
by  the  water.  In  addition,  each  outlet  was  covered  with  two  burlap  bags. 
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the  purpose  being  to  disperse  the  water  Jets  and  to  prevent  the  latter 
from  inducing  wave  motions  in  the  tank  which  would  be  transmitted  back  into 
the  system  as  pressure  fluctuations. 

Instrumentation 

Manometer  System 

In  order  to  facilitate  taking  eighteen  static  pressure  readings  along 
the  test  section,  two  eighteen  tube  manometers  were  constructed  so  that  all 
readings  could  be  obtained  at  one  time.  One  manometer  was  in  a vertical 
position,  and  was  used  when  the  pressure  differences  between  successive  taps 
we re  comparatively  large.  For  low  flow  rates,  when  the  pressure  differences 
were  small,  an  inclined  manometer  was  used.  This  was  positioned  to  give  the 
readings  an  amplification  of  five.  Both  nanometers  were  provided  with  vents 
at  their  bases  so  that  all  air  could  be  removed  from  the  connecting  lines. 
Thermocouples 

Thermocouples  were  used  to  measure  the  entrance  and  exit  water 
temperatures  of  the  test  section,  the  temperature  of  the  tube  wall  at  six 
positions  along  its  length,  and  the  wall  temperature  of  the  small  contraction 
nozzle.  All  of  the  thermocouples  were  made  from  the  same  spool  of  copper- 
constantan  duplex  wire  and  had  welded  junctions.  The  leads  were  brought 
into  an  isothermal  zone  box  and  were  then  connected  through  a selector  switch 
to  a standard  laboratory  potentiometer,  A thermometer  was  suspended  in  the 
storage  tank  as  a check  on  the  average  water  temperature. 

Orifice  Meters 

Standard  ASME  flange-tap  orifice  plates  were  used  to  measure  the  flow 
at  pump  discharge,  and  the  quantity  of  water  removed  through  the  boundary 
layer  suction  slot.  The  flow  in  the  test  section  was  obtained  as  the 
difference  between  these  two  values.  Both  orifice  unions  were  connected  to 
vertical  U-tube  mercury  manometers,  and  by  the  use  of  vents  the  connecting 
lines  were  kept  filled  with  water. 
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Boundary  Layer  Probe 

A boundary  layer  probe,  constructed  as  shown  in  Figure  5,  served  as  a 
second  means  of  locating  the  transition  point.  The  probe  measured  the 
static  pressure  on  the  axis  of  the  tube,  and  the  total  pressure  in  the 
boundary  layer  at  a fixed  distance  from  the  tube  wall  and  in  the  same 
longitudinal  plane  as  the  static  pressure.  The  probe  was  inserted  through 
an  O-ring  seal  at  the  exit  end  of  the  test  pipe,  and  was  calibrated  in 
inches  so  that  it  could  be  positioned  at  any  desired  place  along  the  tube. 


10 


IV,  TEST  PROCEDURE 

At  the  beginning  of  the  day  in  which  the  apparatus  was  being  used,  the 
storage  tank  was  freshly  filled  from  the  city  water  mains,  and  then,  by 
closing  appropriate  valves,  the  entire  flow  was  diverted  through  the  filter 
at  low  velocity  and  circulated  back  into  the  tank.  This  was  continued  for 
almost  an  hour  to  assure  that  essentially  all  of  the  water  had  passed 
through  the  filter  and  was  as  clean  as  possible.  The  main  valves  were  then 
opened  and  a large  flow  allowed  to  circulate  through  the  apparatus  thus 
helping  to  expel  any  air  which  might  be  entrapped.  In  addition  vents  were 
opened  to  remove  air  trapped  in  the  upper  region  of  the  honeycomb  and 
stilling  chamber.  Steam  was  admitted  to  the  jacket  to  heat  the  water  to 
room  temperature  so  that  the  temperature  would  remain  substantially  constant 
during  the  adiabatic  rams.  After  the  water  had  reached  the  proper  condition 
the  steam  was  shut  off  and  the  apparatus  allowed  to  cool  while  the  manometer 
lines  were  carefully  purged  of  entrapped  air.  The  main  flow  and  boundary 
layer  suction  rates  were  then  adjusted  to  the  desired  values  for  the 
experiment. 

Runs  Using  Wall  Static  Pressure  Taps 

The  first  group  of  runs  was  taken  using  the  pressure  readings  of  the 
wall  static  pressure  taps  as  the  means  of  determining  the  transition  point. 
At  low  flow  rates  the  readings  were  taken  from  the  inclined  manometer,  while 
at  higher  flow  rates  when  the  pressure  differences  between  successive  taps 
were  larger,  the  vertical  manometer  was  used.  Immediately  after  taking  an 
adiabatic  run,  the  steam  was  turned  on  without  making  any  changes  in  the 
valve  settings.  Cooling  water  was  admitted  to  the  storage  tank  (which  was 
maintained  at  a constant  level  by  an  overflow  pipe)  to  prevent  the  system 
temperature  from  rising  as  the  result  of  heat  added  through  the  walls  of  the 
test  section.  It  was  found  that  some  control  over  the  size  of  the  tempera- 
ture difference  across  the  wall  could  be  obtained  by  the  amount  of  venting 
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of  the  steam  jacket.  The  more  thoroughly  the  air  was  vented  from  the  jacket 
the  higher  the  temperature  difference  could  be  made,  up  to  about  100°F, 
which  was  the  upper  limit  for  the  low  pressure  steam  employed.  If  high 
temperature  differences  were  used  it  was  found  that  non-condensible  gases 
would  come  out  of  solution  with  the  water,  and  air  bubbles  would  collect  in 
the  manometer  lines,  thereby  resulting  in  erroneous  readings.  When  the 
wall  pressure  taps  were  used  as  a measuring  means,  the  temperature 
differences  were  therefore  kept  small. 

During  the  heated  runs,  cooling  water  was  circulated  through  the  jacket 
of  the  small  contraction  and  regulated  so  that  the  contraction  wall  tempera- 
txire  was  within  a few  degrees  of  the  entering  water  temperature.  In  general 
a warming  up  period  of  at  least  an  hour  was  necessary  during  the  heated  runs 
to  achieve  thermal  equilibrium.  During  this  period  thermocouple  readings 
were  taken  to  aid  in  the  regulation  of  the  cooling  water  flow. 

Runs  Using  the  Boundary  Layer  Probe 

For  the  r*uns  using  the  boundary  layer  probe  it  was  not  necessary  to 
vent  the  manometer  lines,  thus  making  it  possible  to  reduce  the  starting 
time  of  the  equipment.  The  connections  of  the  probe  to  the  inclined 
manometer  were  opened,  the  connections  to  the  static  wall  taps  were  closed, 
and  the  lines  running  from  the  probe  were  thoroughly  purged  of  air.  The 
first  reading  was  usually  taken  one  inch  downstream  fro'--:  the  boundary  layer 
suction  slot,  and  successive  readings  were  obtained  at  various  intervals 
along  the  pipe  by  sliding  the  probe  down  the  tube. 

For  a particular  main  stream  flow,  runs  were  made  with  various  rates 
of  boundary  laver  suction  to  determine  whether  or  not  the  location  of  the 
transition  point  was  sensitive  to  the  suction  rate.  The  runs  with  heat 
applied  werejthen  taken  Within  a range  whefe  transition  was  insensitive  to 
the  suction  rate,  so  that  any  change  observed  could  not  be  attributed  to 
slight  differences  in  suction  rate  between  the  adiabatic  runs  and  the 
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heated  runs.  Since,  in  the  case  of  the  boundary  layer  probe,  all  of  the 
connecting  lines  ran  down  the  interior  of  the  test  pipe  and  were  therefore 
kept  cool  by  the  water  flow,  the  difficulty  of  non-condensible  gases 
collecting  in  the  lines  was  not  present  and  higher  temperature  differentials 
could  be  applied. 

At  the  end  of  the  running  period  the  apparatus  was  completely  drained 
to  help  minimize  corrosion  or  electrolytic  action  between  dissimilar  metals 


within  the  system 
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Vj_  EXPERIMENTAL  RESULTS 

Methods  Used  to  Determine  the  Location  of  Transition 
Static  Pressure  Taps 

When  a fluid  flows  through  a pipe  entry,  a laminar  boundary  layer  grows 
which  eventually  enters  a transition  zone  and  becomes  turbulent.  The  growth 
of  the  boundary  layer  is  accompanied  by  momentum  changes  caused  by  altera- 
tions in  the  velocity  profile  as  the  flow  pattern  develops.  Thus  the 
difference  in  static  pressure  between  two  points  along  the  pipe  is  a 
measure  not  only  of  the  wall  friction  drop,  but  also  includes  the  change 
due  to  momentum  effects.  The  entire  static  pressure  change  between  two 
points  is  expressed  conveniently  in  dimensionless  form  as  a local  apparent 
friction  factor,  defined  as  follows: 

4fAPP  E dp/d(x/D) 

if>v2 

In  the  laminar  zone  the  local  apparent  friction  factor  varies  inversely 
as  the  one-half  power  of  the  length  Reynolds  number,  and  thus  decreases 
along  the  length  of  the  pipe.  During  transition  the  friction  factor  rises 
quite  rapidly  (it  may  drop  at  first  under  some  circumstances)  as  the 
boundary  layer  thickens,  and  then  begins  to  decrease  slowly  after  the 
turbulent  layer  has  formed.  By  considering  a plot  of  4fApp  vs.  Rex  oq, 
logarithmic  coordinates,  the  transition  point  can  be  determined  where  the 
curve  begins  to  deviate  from  its  original  slope  of  negative  one-half. 
Boundary  Layer  Probe 

The  boundary  layer  probe  was  constructed  to  measure  the  static  oressure 
along  the  axis  of  the  tube  in  the  core  of  the  flow,  as  well  as  the  total 
pressure  in  the  boundary  layer  at  a fixed  distance  from  the  wall.  The 
difference  between  the  two  readings  obtained  is  a direct  indication  of  the 
local  velocity  at  a fixed  distance  from  the  wall.  Owing  to  the  severe 
change  in  shape  of  the  velocity  profile  at  transition,  this  technique 
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provided  a convenient  means  for  locating  the  transition  point. 

As  the  boundary  layer  thickness  increases  along  the  tube  in  the 
laminar  zone,  the  probe,  located  at  increasing  axial  lengths  but  at  a 
constant  distance  from  the  wall,  becomes  positioned  proportionately  deeper 
within  the  layer.  Thus  it  is  placed  in  a region  of  lower  velocity,  and  a 
lower  dynamic  pressure  is  recorded  as  the  difference  between  total  and  static 
pressure  readings.  In  the  laminar  region,  the  curves  of  velocity  pressure 
versus  length  Reynolds  number  have  a negative  slope  similar  to  the  local 
friction  factor  curves  previouslv  discussed.  To  further  demonstrate  this 
an  approximate  velocity  profile  for  laminar  flow  over  a flat  plate  was 
selected  from  reference  18,  and  calculations  performed  as  indicated  in 
Appendix  B.  This  resulted  in  the  illustrative  curve  shown  in  Figure  6. 

The  calculated  curve  for  the  laminar  zone  shows  the  general  trend  of 
negative  slope  and  agrees  fairly  well  with  the  shape  obtained 
experimentally  in  this  region. 

When  the  boundary  layer  becomes  turbulent  the  resulting  change  in  the 
shape  of  the  velocity  profile  is  reflected  as  an  increase  in  the  velocity 
pressure  determined  by  the  probe.  The  transition  point  is  thus  located  as 
a well  defined  minimum  on  the  curve  of  velocity  pressure  versus  length 
Reynolds  number. 

Since  the  boundary  layer  probe  was  moved  along  the  tube  to  obtain 
readings  at  various  points,  all  of  the  readings  were  taken  bv  the  same 
instrument.  Thus  any  small  errors  inherent  in  the  instrument,  (e.g.,  owing 
to  irregularities  in  the  pressure  tap  openings),  remained  substantially 
constant  for  all  readings  and  did  not  cause  any  scatter  of  the  data  around 
the  curves.  Since  the  velocities  used  were  fairly  sizeable,  the  velocity 
pressure  readings  obtained  were  much  larger  than  those  due  to  friction 
drop.  This  allowed  a wider  margin  of  experimental  reading  error  without 
causing  noticeable  scatter  of  the  data.  For  these  reasons  the  results 
obtained  with  the  boundary  layer  probe  produced,  in  general,  better 
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experimental  curves  than  the  wall  pressure  tap  method. 

One  question  which  immediately  arises  is  whether  the  disturbance  of 
the  boundary  layer  probe  in  the  test  section  may  be  transmitted  upstream 
and  possibly  cause  transition  to  occur  earlier  than  if  the  probe  were  not 
present.  To  studv  this  effect  a run  was  made  in  which  transition  was 
determined  by  both  methods,  the  wall  static  pressure  readings  being  taken 
after  the  probe  had  been  withdrawn  beyond  the  end  of  the  test  section. 

The  two  curves,  shown  in  Figures  7 and  9,  both  indicate  transition  at  the 

5 

same  length  Reynolds  number  of  about  4.7  x 10  . This  demonstrates  that  the 
introduction  of  the  probe  did  not  significantly  affect  the  position  of  the 
transition  point. 

Graphical  Presentation  of  the  Data 

Either  the  local  apparent  friction  factor  or  the  velocity  pressure  in 
the  boundary  layer  has  been  plotted  against  the  length  Reynolds  number  on 
logarithmic  coordinates.  The  value  of  the  average  friction  factor  measured 
between  two  taps  has  to  be  plotted  at  a length  Reynolds  number  corresponding 
to  some  mean  value  of  x between  the  two  points.  In  the  laminar  zone,  as 
shown  in  reference  6,  the  appropriate  mean  value  is  the  average  of  the 
arithmetic  and  geometric  means  of  the  axial  distances  of  the  two  taps  from 
the  boundary  layer  suction  slot  at  the  entrance  of  the  test  section. 

Two  different  length  Reynolds  numbers  have  been  used  for  plotting  the 
curves.  The  free  stream  Reynolds  number,  Rex,  was  based  on  the  average 
main  stream  water  temperature  in  the  core  of  the  flow.  The  mean  boundary 
layer  Reynolds  number,  Re^,  was  based  on  the  mean  boundary  layer  tempera- 
ture; the  latter  was  taken  as  the  average  between  the  wall  and  main  stream 
temperatures. 

With  the  exception  of  Figure  9,  the  data  placed  on  each  figure  were 
all  obtained  during  the  same  day.  This  was  to  assure  that  the  runs  shown 
on  each  graph  could  be  directly  compared. 


16 


Discussion  of  the  Data 

Frlotlon  Factor  Data 

Effect  of  Dirt  Particles  on  Transition:  Figures  9 to  12  illustrate 

the  performance  of  the  system  for  a variety  of  diameter  Reynolds  numbers 
ranging  from  115,000  to  167,000,  with  transition  determined  from  measure- 
ments of  local  apparent  friction  factor.  A predominant  characteristic  which 
the  curves  display  is  that,  as  successive  runs  were  taken,  the  value  of  the 
transition  Reynolds  number  decreased  from  a maximum  of  300,000  in  Figure  9 
to  400,000  in  Figure  11.  This  is  especially  evident  from  a comparison  of 
Figures  10  and  11  which  have  essentially  the  same  diameter  Reynolds  numbers 
and  percent  boundary  layer  suction,  but  different  transition  Reynolds 
numbers.  This  characteristic  of  decreasing  transition  Revnolds  number  on 
successive  days  existed  throughout  the  entire  period  of  experimentation, 
and  was  attributed  to  the  collection  of  fine  particles  of  foreign  matter  on 
the  surfaces  of  the  inlet  nozzles  and  test  section.  This  increased  the 
disturbance  level  of  the  system,  thereby  causing  transition  to  occur  earlier. 
This  conclusion  was  verified  by  the  fact  that  the  transition  point  could  be 
restored  to  its  earlier  position  by  giving  the  inlet  sections  and  test  pipe 
a thorough  cleaning. 

Effect  of  Heating  on  Transition:  Figures  10,  11,  and  12,  each  contain 

two  curves,  one  for  adiabatic  conditions,  and  the  other  with  steam  supplied 
to  produce  a small  temperature  difference  across  the  boundary  laver.  The 
curves  show  that  heating  caused  no  appreciable  change  in  the  free  stream 
transition  Reynolds  number,  but  that  the  mean  boundary  layer  transition 
Revnolds  number  was  increased  slightly. 

Data  From  Boundary  Layer  Probe 

Figure  13  shows  the  first  data  taken  with  the  boundary  layer  probe 
immediately  after  the  system  had  been  given  a thorough  cleaning.  It  is 
noticed  that  the  decrease  in  disturbance  level  due  to  cleaning  has  doubled 
the  transition  Reynolds  number  as  compared  to  the  previous  run. 
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Effect  of  Boundary  Layer  Suction  on  Transition:  Figures  14,  17,  18, 

19,  and  20  show  experimental  results  for  a range  of  diameter  Reynolds 
numbers  from  30,000  to  120,000,  For  each  adiabatic  set  of  curves  with 
constant  Rep  there  are  shown  various  rates  of  boundary  layer  suction  flow. 
These  curves  illustrate  that  beyond  a certain  percentage  the  suction  rate 
did  not  have  an  appreciable  effect  on  the  location  of  transition. 

When  the  boundary  layer  suction  rate  was  zero,  on  the  other  hand,  the 
transition  point  was  displaced  considerably  upstream.  This  is  partly 
explained  by  the  fact  that  the  boundary  layer  formed  in  the  inlet  section 
was  not  being  removed,  thereby  producing  the  effect  of  extending  the  length 
of  the  pipe  upstream,  and  partly  by  the  fact  that  the  boundary  layer 
suction  slot  itself  was  causing  a disturbance  under  these  circumstances. 

Effect  of  Heating  and  Accumulation  of  Foreign  Particles  on 
Transition:  Compared  with  the  adiabatic  curves  in  Figures  14,  17,  18,  19, 

and  20,  the  heating  curves  exhibit  a slightly  smaller  free  stream  Reynolds 

y 

number  of  transition,  but  a slightl7r  greater  mean  boundary  layer  Reynolds 
number  of  transition.  The  heating  curves  are  also  positioned  above  the 
adiabatic  curves,  thus  giving  evidence  that  the  velocity  profile  near  the 
wall  has  been  modified  by  heat  transfer. 

Throughout  this  series  of  runs  it  may  again  be  noted  that  there  was  a 
continuous  decrease  in  the  transition  Reynolds  number  due  to  the  accumula- 
tion of  dirt  and  foreign  matter  in  the  inlet  sections.  During  the 
experimental  work  the  location  of  the  transition  point  remained  constant 
during  any  particular  day.  An  example  of  this  is  shown  in  Figure  25  where 
runs  68  and  72  show  transition  at  the  same  location  although  one  run  was 
taken  in  the  morning  and  the  other  late  in  the  afternoon.  The  transition 
Rex  appeared  to  decrease  from  day  to  day  and  especially  over  week-ends 
when  the  drained  system  was  allowed  to  stand  for  a few  days.  This  suggests 
that  the  foreign  matter  was  probably  salts  deposited  by  evaporation  of 
water  left  within  the  system.  The  powdery  white  deposits  found  on  the 
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interior  surfaces  of  the  system  when  it  was  cleaned  seemed  to  he  of  this 
nature . 

Since  it  was  desirable  to  eliminate  as  much  as  possible  the  effect  of 
stream  disturbances,  the  system  was  given  another  cleaning  before  obtaining 
run  60  in  Figure  21.  This  caused  the  transition  Reynolds  number  to  increase 
to  twice  the  value  indicated  in  Figure  20.  Figures  21  and  22  again 
demonstrate  that  heating  was  ineffective  in  increasing  the  free  stream 
length  Reynolds  number  of  transition,  but  that  the  mean  boundary  layer 
Reynolds  number  of  transition  was  slightly  increased. 

Because  the  cleaning  had  only  increased  the  adiabatic  transition 
Reynolds  number  to  about  500,000,  it  was  decided  to  give  the  system  still 
another  cleaning,  as  higher  values  had  been  reached  previously.  Before 
run  64  in  Figure  23  was  obtained  the  system  had  been  given  a very  thorough 
cleaning  and  the  nozzles,  stilling  chamber,  and  test  section  were  brightly 
polished.  The  joint  between  the  two  contractions  was  further  honed  to 
assure  its  being  especially  smooth. 

Figures  23  and  24  show  the  first  group  of  runs  taken  after  cleaning, 
and  it  is  seen  that  the  adiabatic  transition  Reynolds  number  was  restored 
to  850,000.  It  is  also  evident  that  heating  produced  a rather  pronounced 
decrease  in  the  free  stream  transition  Reynolds  number,  but  a slight 
increase  in  the  mean  boundary  layer  transition  Reynolds  number.  Figures 
25,  26,  27,  and  28  show  the  effects  of  heating  with  different  temperature 
gradients  across  the  boundary  layer.  In  all  cases  heating  appeared  to 
cause  a small  decrease  in  the  free  stream  transition  Reynolds  number  but  an 
increase  in  the  mean  boundary  layer  transition  Reynolds  number.  The 
increase  in  the  latter  was  generally  more  pronounced  as  the  temperature 
difference  was  increased. 

The  Addition  of  Detergent  to  the  System:  During  a previous  experimental 

work  in  which  lucite  contraction  nozzles  were  employed  in  a water  flow 
system,  it  was  observed  that  tiny  air  bubbles  would  adhere  to  the  interior 
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surfaces  of  the  nozzles  and  cause  large  reductions  in  the  transition 
Reynolds  number  due  to  their  effective  roughness.  In  the  present  apparatus 
it  was  thought  that  a small  accumulation  of  grease  on  the  nozzle  surfaces 
would  facilitate  the  adherence  of  such  bubbles  and  might  possibly  account 
for  the  observed  day  to  dav  decrease  of  transition  Revnolds  number. 
Detergent  was  therefore  added  to  the  water  during  run  63  to  eliminate 
any  grease  film  and  bubble  collection,  and  in  addition  to  help  prevent 
the  accumulation  of  dirt  particles.  Finrure  29  shown  that  the  addition  of 
the  detergent  was  ineffective  in  displacing  the  position  of  transition* 

This  result  does  not  eliminate  the  possibility  of  air  bubbles  coming 
out  of  solution  and  remaining  in  suspension,  with  the  net  effect  of  causing 
disturbances  which  might  lead  to  premature  transition. 
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VI.  SIGNIFICANCE  OF  RESULTS 

The  results  of  these  experiments  have  clearly  demonstrated  that  in  the 
case  of  internal  water  flows,  the  length  Reynolds  number  of  transition  is 
very  strongly  dependent  upon  the  stream  disturbance  level.  Sven  very  small 
deposits  of  foreign  matter  along  the  walls  were  sufficient  to  reduce  the 
value  of  Rex  at  transition  by  fifty  percent.  Sven  with  special  precautions 
it  was  difficult  to  attain  values  of  Rex  at  transition  as  high  as  850,000, 
and  in  the  usual  water  system  of  a practical  nature  it  would  be  difficult 
to  achieve  such  a low  disturbance  level. 

For  values  of  disturbance  level  present  in  this  work,  corresponding 
to  values  of  Rex  at  transition  ranging  from  850,000  to  250,000,  no  increase 
in  the  free  stream  length  Revnolds  number  of  transition  was  observed  when 
heat  was  applied,  and  in  fact  a slight  decrease  occurred  in  some  instances. 
However,  heating  consistently  caused  an  increase  in  the  mean  boundary  layer 
length  Reynolds  number  of  transition.  Since  in  design  work  it  is  the  free 
stream  Reynolds  number  which  is  physically  significant  in  predicting  the 
location  of  transition,  it  can  be  said  that  with  normal  disturbance  levels 
present  in  practical  situations,  heating  would  be  of  no  significance  in 
extending  the  laminar  regime  of  flow. 

These  experiments  have  not  conlusively  proven  that  heating  is 
ineffective  in  delaying  transition,  for  if  the  stream  disturbance  level 
could  be  sufficiently  reduced,  then  perhaps  a heatin'?  effect  would  be 
realized.  However,  as  shown  by  the  experimental  results,  it  is  extremely 
difficult  to  obtain  streams  of  sufficiently  low  turbulence  level,  and  so 
it  is  doubtful  that  heating  could  be  utilized  as  a practical  means  for 

delaying  transition  except  perhaps  under  unusual  circumstances, 

i * 1 > 

At  present  the  mathematical  theory  of  stability  has  in  fact  not  been 

developed  for  the  case  of  liquid  flows  with  heating  or  cooling,  so  we  can 
only  speculate  as  to  what  the  results  may  be  by  drawing  an  analogy  with  the 
results  for  air  flows.  Since  for  a liquid  the  behavior  of  viscosity  as  a 
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function  of  temperature  Is  opposite  to  that  of  a gas,  it  would  be  expected 
that  the  Influence  of  temoerature  on  the  shape  of  the  velocity  profile  in 
a liquid  flow  would  also  be  opposite  to  that  of  a gas  flow.  Thus,  while 
heating  tends  to  destabilize  a gas  flow,  we  mi<?ht  expect  it  to  cause  a 
liquid  flow  to  become  more  stable.  Fowever,  this  analogy  could  only  apply 
for  predicting  the  behavior  of  the  instability  point  in  a stream  containing 
vanishingly  small  disturbances. 

Assuming  that  the  analogy  is  valid  for  liquid  flows,  the  present  work 
provides  two  results:  1)  In  a practical  case  where  the  disturbances  are 

not  vanishingly  small,  their  effect  may  be  very  significant  and  perhaps 
actually  govern  the  process.  2)  The  results  provide  additional  evidence 
that  in  many  circumstances  it  is  not  possible  to  extend  predictions  for  the 
location  of  the  instability  point  to  that  of  the  transition  point. 

It  must  be  noted  that  there  are  really  two  effects  entering  when  heat 
is  acplied  to  the  wall  of  the  test  section.  One  is  the  change  in  shape  of 
the  velocity  profile  resulting  from  the  viscosity  gradient  produced  by 
heating  the  boundary  layer.  This  may  cause  an  increase  or  decrease  in 
stability  which  may  in  turn  change  the  position  of  transition.  The  other 
is  merely  the  change  in  the  mean  viscosity  of  the  boundary  layer.  When  a 
liquid  boundary  layer  is  heated,  the  mean  viscosity  decreases,  and  hence 
the  effective  Reynolds  number  in  the  layer  is  increased  for  the  same  flow 
velocity.  This  per  se  will  tend  to  cause  transition  to  occur  earlier. 
Evidently,  from  the  experimental  results,  if  stability  is  increased  by 
heating,  the  viscosity  effect  is  sufficiently  large  to  compensate  and  cause 
the  net  result  to  be  essentially  zero  or  slightly  negative.  Actually,  in 
this  respect,  the  mean  boundary  layer  Reynolds  number  is  a better  indica- 
tion of  what  occurs  within  the  boundarv  laver  itself.  TJhe  delay  in  transi- 
* 

tion  in  relation  to  that  parameter  indicates  that  an  increase  in  stability 
was  present,  although  not  sufficiently  large  to  be  of  practical  importance. 

After  the  mathematical  theory  has  been  investigated  for  liquids  to  an 
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extent  comparable  to  the  present  knowledge  for  air,  two  additional  things 
remain  to  be  determined  before  comparisons  with  experiments  can  be  made 
with  a better  degree  of  certainty#  First  the  effect  of  large  disturbances 
should  be  dealt  with,  that  is,  the  non-linear  effects  of  the  governing 
equations  must  be  considered.  Secondly,  the  zone  of  unstable  wave  amplifi- 
cation which  occupies  the  interval  between  the  instability  point  and  the 
transition  point  should  be  investigated,  including  the  effect  of  tempera- 
ture gradients  on  amplification.  At  least  a qualitative  answer  to  the 
latter  question  is  necessary  before  instability  criteria  can  be  correlated 
with  the  behavior  of  transition. 
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VII.  RECOMMENDATIONS  FOR  FUTURE  'YORK 

If  further  work  is  performed  with  this  type  of  system  an  improvement 
would  result  if  the  system  were  completely  closed  and  contained  distilled 
water.  This  would  help  prevent  dirt  and  insoluble  salts  from  collecting  on 
the  walls  of  the  system,  and  in  addition  would  eliminate  the  possible 
appearance  of  air  bubbles.  Since  cooling  water  would  not  be  admitted 
during  the  heated  runs  under  these  circumstances,  refrigeration  would  be 
necessary  to  keep  the  system  temperature  from  rising. 

Since  the  boundary  layer  probe  proved  successful  as  a means  for 
determining  transition,  the  use  of  wall  pressure  taps  could  be  eliminated. 
This  would  prevent  any  possible  disturbance  that  the  taps  themselves  may 
induce  and  would  also  simplify  the  apparatus  construction. 

It  would  be  worthwhile  to  investigate  the  behavior  of  the  boundary 
layer  under  the  effect  of  larger  temperature  gradients.  This  could  be 
accomplished  by  either  pressurizing  the  steam  jacket  to  accommodate  high 
temperature  steam,  or  by  maintaining  a lower  water  temoerature  by  the 
use  of  refrigeration. 

It  would  be  of  interest  to  try  cooling  the  liquid  boundary  layer.  In 
this  instance  the  increase  of  viscosity  in  the  boundary  layer  produces  the 
effect  of  decreasing  the  Reynolds  number  in  that  region,  and  this  would 
tend  to  move  transition  downstream.  On  the  other  hand,  the  boundary  layer 
velocity  profile  would  tend  to  develop  an  inflection,  thus  becoming  more 
unstable. 
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VIII.  FRICTION  FACTOR  IN  THE  LAMINAR 
flNTtiY  OF  A SMOCtTTpOT 

Brief  Summary  of  Previous  '.York 

Various  theories  predicting  the  behavior  of  the  boundary  laver  for 
adiabatic,  incompressible  flow  in  the  laminar  entry  of  a tube  have  been 
formulated  by  Boussinesq^  ^ , Schiller^®),  and  Langhaar^ ^ . Atkinson  and 
Goldstein^)  have  solved  the  differential  equations  by  use  of  series 
expansions  and  found  expressions  valid  near  the  entrance  of  the  pipe.  This 
solution  was  then  joined  to  a second  exoression,  valid  at  larger  distances 
from  the  entrance,  derived  from  an  extension  of  the  method  used  by 
Boussinesq.  A simple  theory,  using  momentum  integral  methods,  has  been 
derived  bv  Kline  and  Shapiro^6)  for  the  region  where  the  boundary  layer  is 
very  thin.  They  also  present  a set  of  curves  comparing  the  calculated 
results  for  all  these  methods  with  their  experimental  data. 

Solution  for  Region  Where  S/R  is  Very  Small 

The  development  of  the  boundary  layer  in  a tube  is  complicated  by  the 
fact  that  its  formation  affects  the  velocity  in  the  frictionless  core  of  the 
flow.  The  resulting  acceleration  of  the  core  produces  a pressure  gradient 
along  the  pipe  which  in  turn  affects  the  boundary  layer  formation.  The 
basic  physical  nature  of  the  process  suggests  that  a method  of  successive 
approximations  might  be  used  to  obtain  a solution.  In  this  scheme  the 
continuity  equation  is  applied  to  obtain  an  approximation  to  the  core 
velocity  along  the  tube.  The  latter  is  then  substituted  in  the  approximate 
boundary  layer  calculation  method  of  Thwaites^^,  and  of  Rott  and 
Crabtree^l^),  to  obtain  the  displacement  thickness  as  a function  of  distance 
along  the  tube.  This  displacement  thickness  is  then  substituted  back  into 
the  continuity  relation  to  obtain  a better  approximation  to  the  core  velocity, 
and  the  whole  operation  is  repeated. 

The  calculations,  as  carried  out  in  Appendix  C,  result  in  the 


expression, 
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^>fpl  . 4-  III  ft,.-...' 

l Be,  ft,'  I 

which  is  valid  in  the  range  where  Re^/Rei)2  is  less  than  10"^,  within  which 
good  convergence  is  obtained.  For  Re^Rep2  equal  to  10“3,  S/R  is 
approximately  0.3, 

-3 

For  values  larger  than  10  the  boundary  layer  has  become  so  thick  that 
three  dimensional  effects  enter  and  the  method  breaks  down.  Also  since  the 
solution  has  not  been  constrained  to  converge  eventually  to  a fully- 
developed  laminar  profile  it  would  be  expected  to  be  valid  only  in  the  entry 
region  of  the  pipe  where  the  boundary  layer  is  still  relatively  thin 
compared  to  the  pipe  diameter.  A comparison  of  the  computed  values  with 
experimental  data  is  shown  in  Figure  30. 

Solution  for  Entire  Laminar  Entry  Region 
A second  solution  for  the  apparent  friction  factor  in  a laminar  tube 
entry  was  obtained  by  substituting  an  approximate  expression  for  the  velocity 
profile  into  the  eqtiations  of  continuity  and  momentum.  The  profile  chosen 
was  in  the  form  u/U  * Function  £(y/S),  ( &/R )j  , and  was  constrained  to 
approach  a parabolic  shape  as  the  boundary  layer  thickness  became  equal 
to  the  tube  radius.  This  method  yielded  a solution  which  applies  throughout 
the  whole  entry  region,  and  is  in  close  agreement  with  both  the  experimental 
data  of  Kline  and  Shapiro^®)  and  with  the  expression  stated  in  the  preceding 
paragraph.  The  computations  are  fully  explained  in  Appendix  D,  and  the 
results  are  shown  in  Figures  31  to  34. 
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IX.  CONCLUSIONS 

1.  For  disturbance  levels  producing  length  Reynolds  numbers  of  transition 
up  to  850,000,  the  extent  of  the  laminar  regime  for  water  flowing  in  a pipe 
entry  was  found  to  be  either  unaffected  or  slightly  decreased  when  the  pipe 
walls  were  heated, 

2.  The  length  Reynolds  number  of  transition  for  water  flows  was  found  to 
be  very  sensitive  to  small  disturbances  such  as  might  be  caused  by  small 
roughnesses  on  the  walls. 

3.  Heating  caused  small  increases  in  the  mean  boundary  layer  length 
Reynolds  number  of  transition,  which  indicates  an  increased  stability. 
However,  the  increases  produced  were  not  sufficiently  large  to  be  of 
practical  importance, 

4.  There  is  need  for  further  mathematical  study  of  the  unstable  wave 
amplification  which  follows  the  instability  point  and  induces  turbulence, 
before  a clear  understanding  of  the  relation  between  the  instability  and 
transition  points  can  be  obtained. 


27 


REFERENCES 


1.  Boussinesq,  J.,  "Hydrodynamique, " Comptes  Rendus,  Vol,  110,  1890, 

pp.  1160-1238;  Vol,  113,  1891,  pp.  9-49. 

2.  Dryden,  H.  L.  and  Abbot,  I,  H, , "The  Design  of  Low-Turbulence  Wind 

Tunnels,"  NACA  T.  N.  No.  1755.  Nov.  1949. 

3.  Dryden,  H.  L.  and  Schubauer,  G.  B.,  "The  Use  of  Damping  Screens  for  the 

Reduction  of  Wind-Tunnel  Turbulence,"  Journal  of  the  Aero.  Scl.. 

Vol,  14,  April  1947,  pp.  221-228. 

4.  Eckert,  B.  and  Pfluger,  F.,  "The  Resistance  Coefficient  of  Commercial 

Round  Wire  Grids,"  NACA  T.  WL  No^  1003.  Jan.  1942. 

5.  Goldstein,  S.,  Modern  Developments  in  Fluid  Dynamics.  Oxford  University 

Press,  1938. 

6.  Kline,  S.  J.  and  Shaoiro,  A.  H.,  "The  Effect  of  Cooling  on  Boundary 

Layer  Transition  in  a Gas,"  Final  Report  on  DIC  Froject  3-6927 
at  M.I.T.  Submitted  to  NACA  under  Contract  NAw-6125,  Sept,  30,  1952. 

7.  Langhaar,  H.,  "Steady  Flow  in  the  Transition  Length  of  a Straight  Tube," 

Journal  of  Applied  Mech.,  Vol,  9,  No.  2,  June  1942,  pp.  A-55  to  A-58. 

8.  Lees,  L.,  "The  Stability  of  the  Laminar  Boundary  Layer  in  a Compressible 

Fluid,"  NACA  Report  No.  876.  1947. 

9.  Lees,  L.  and  Lin,  C,  C.,  "Investigation  of  the  Stability  of  the  Laminar 

Boundary  Layer  in  a Compressible  Fluid, " NACA  T.  N,  No.  1115. 

Sept.  1946. 

10.  Liepmann,  H.  W.  and  Fila,  G.,  "Investigations  of  Effects  of  Surface 

Temperature  and  Single  Roughness  Elements  on  Boundary  Layer 
Transition,"  NACA  T\  No^  1196,  April  1947. 

11.  Lin,  C.  C.,  "On  the  Stability  of  Two-Dimensional  Parallel  Flows," 

Quarterly  of  App . Math. , Part  I,  Vol.  Ill,  No.  2,  July  1945, 
pp.  117-142;  Part  I 1,  Vol.  Ill,  No.  3,  Oct.  1945,  pp.  219-234; 

Part  III,  Vol.  Ill,  No.  4,  Jan.  1946,  pp.  277-301. 

12.  McCullough,  G.  B.  and  Gault,  D.  E.,  "An  Experimental  Investigation  of 

the  NACA  63i~012  Airfoil  Section  with  Leading-Edge  and  Midchord 
Suction  Slots,"  NACA  T^  N^  2041,  Feb.  1950. 

13.  Rott,  N.  and  Crabtree,  L.  F.,  "Simplified  Laminar  Boundary-Layer 

Calculations  for  Bodies  of  Revolution  and  for  Yawed  Wings, " 

Journal  of  the  Aero  Scl. . Vol,  19,  No.  8,  Aug.  1952. 

14.  Rouse,  H.  and  Hassan,  M.  M.,  "Cavitation-Free  Inlets  and  Contractions," 

Mechanical  Engineering.  Vol.  71,  No.  3,  March  1949,  pp.  213-216. 

15.  Scherrer,  R.,  "Boundary-Layer  Transition  on  a Cooled  20°  Cone  at  Mach 

Numbers  of  1.5  and  2.0,"  NACA  T.  N . No.  2131.  July  1950. 


28 


16.  Schiller,  L.,  "Die  Entwicklung  der  Laminaren  Geschwindigkeitverteilung 

und  ihre  Bedeutung  fur  Zahigkeitsmessungen,  " Zeitschrlft  fur 
Angewandte  Mathematic  und  Me chan lk,  Vol.  2,  1922,  p,  96. 

17.  Schlichting,  H.,  "Amplltudenverteilung  und  Enerbilanz  der  kleinen 

Storungen  bei  der  Plattenstromung, " Nachr.  d.  Ges.  cU  Wls3_.  zu 
Gottingen,  Math.  Phys.  Kl,,  Neue  Folge,  Bd.  1,  Nr.  4,  1935, 
pp . 47-43. 

13.  Schlichting,  H.,  Grenzschicht-Theorie , Karlsruhe:  Verlag  G.  Braun, 

1951. 

19.  Schubauer,  G.  B.  and  Skramstad,  H.  K.,  "Laminar-Boundary  Layer 

Oscillations  and  Transition  on  a Flat  Plate, " NACA  Report 
No,  909,  1943. 

20.  Schubauer,  G.  B.,  Spangenberg,  W.  G.,  and  Kebanoff,  P.  S.,  "Aerodynamic 

Characteristics  of  Damping  Screens,"  NACA  T.  N.  No . 2001,  Jan. 
1950. 

21.  Thwaites,  B.,  "Approximate  Calculation  of  the  Laminar  Boundary  Layer," 

The  Aeronautical  Quarterly  of  the  Royal  Aero.  Soc. , Vol.  1, 

Nov.  1949,  pp.  245-280. 

22.  Tollmien,  W.,  "Uber  die  Entstehung  der  Turbulenz, " Nachr.  d.  Ges.  d^ 

Wiss.  zu  Gottingen,  Math-Physik  Kl.,  1929,  pp.  21-44 .( translated 
as  NACA  M^  No^  609.  1931.) 

23.  Tollmien,  W.,  "Ein  allgemeines  Kriterium  der  Instabilitat  Laminarer 

Geschwindigkeitverteilungen,  " Nachr.  d.  Ges , d^  'Visa . zu 
Gottingen,  Math-Physik  Kl.,  Neue  Folge,  Bd.  1,  Nr.  5,  1935, 
pp.  79-144.  (translated  as  NACA  T.  M.  No.  792,  1936.) 

24.  Van  Dreist,  E.  R.,  "Calculation  of  the  Stability  of  the  Laminar 

Boundary  Layer  in  a Compressible  Fluid  on  a Flat  Plate  with 
Heat  Transfer,"  Journal  of  the  Aero.  Scl. , Vol.  19,  No.  12, 

Dec.  1952,  pp.  801 -812. 

25.  "Fluid  Meters,  Their  Theory  and  Application,"  ASME  Research 

Publication,  Fourth  Edition,  1937. 


I 


29 


LIST  OP  SYMBOLS 


A 

D 

4fAPP 


47APP 


H 

m 


cross  sectional  area  of  the  test  tube 
diameter  of  the  test  tube 

local  apparent  friction  factor  as  defined  by  the  equations 

4fAPP  i dp/d(x/D) 


mean  apparent  friction  factor  up  to  a section  x = L,  defined 
by  the  equation: 


the  ratio  of  V*  to  0 

a parameter  in  Thwaites'  method  defined  by: 


^STATIC 

P 


m.-  -UV  u'.iU 

^ lx 

static  pressure  in  the  test  section  at  distance  x from 
the  boundary  layer  suction  slot 


PtqtaL  total  pressure  at  x in  the  boundary  layer  at  a fixed  small 
^ distance  from  the  tube  wall 


Q 


volume  rate  of  flow  through  the  test  section 


R pipe  radius 

r radial  coordinate  in  the  test  section  measured  from  the 

tube  axis 


ReD 


R0x 


R®Xm 


To 


diameter  Reynolds  number  based  on  the  diameter  of  the  test 
section  and  on  the  free  stream  fluid  properties,  P*W/u. 


length  Reynolds  number  based  on  the  length  from  the  boundary 
layer  suction  slot  and  on  the  free  stream  properties. 


length  Reynolds  number  based  on  the  length  from  the  boundary 
layer  suction  slot  and  on  fluid  properties  taken  at  the 
average  of  the  free  stream  and  wall  temperatures, 

free  stream  water  temperature  in  the  test  section 


Tw  wall  temperature  of  the  test  section 

U velocity  in  the  frictionless  core  of  the  pipe  flow 
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U„ 

u 

V 

x 

J 

& 

e 

x 


main  stream  velocity  for  flow  over  a flat  plate 

axial  velocity  at  any  location  x and  y 

mean  flow  velocity  through  the  test  section,  Q/k 

axial  distance  in  the  test  section  measured  from  the  boundary 
layer  suction  slot 

radial  coordinate  in  the  test  section  measured  from  the  wall 
a constant 
a constant 

boundary  layer  thickness 

displacement  thickness  of  the  boundary  layer 
momentum  thickness  of  the  boundary  layer 
Pohlhausen  form  parameter  defined  by  the  equation: 

x » <ip 

A* 

coefficient  of  viscosity 
kinematic  viscosity, 
mass  density 

dimensionless  length  parameter,  Re^/Rej^ 
shear  stress  at  the  wall 


I i 


31 


TABLE  I 

LOCATION  OP  STATIC  PRESSURE  TAPS 
AEONS  TEST  SECTION 


Tap  Number 


1 

2 

3 


Distance  Prom 
Suction  Slot 

inches 

1 1/4 

3 

5 


4 

5 

6 


7 

9 

11 


7 


13 


8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 


15 

18 

21 

25 

29 

34 

40 

48 

57 

68 

80 


t 


t 


TABLE  II 


THERMOCOUPLE  LOCATIONS 


Thermocouple  Number 
1 
2 

3 

4 

5 
5 

7 

8 
9 

10 


Location 

Interior  of  honeycomb. 

Boundary  layer  suction 
collection  header. 

Wall  of  the  small 
contraction  nozzle. 

Wall  of  test  section, 
2"  from  suction  slot. 

Wall  of  test  section, 
8"  from  suction  slot. 

Wall  of  test  section, 
14"  from  suction  slot. 

Wall  of  test  section, 
24"  from  suction  slot. 

Wall  of  test  section, 
41"  from  suction  slot. 

Wall  of  test  section, 
70"  from  suction  slot. 

Test  section  exit. 
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APPENDIX  A 
Analysis  of  Data 

The  flow  in  the  main  pipe  line  leading  to  the  test  section  and 
boundary  layer  suction  slot,  and  the  flow  in  the  boundary  layer  suction 
were  both  calculated  from  the  pressure  drop  across  standard  sharp  edge 
orifice  plates.  The  orifice  coefficients  used  were  obtained  from  the 
tables  of  values  supplied  by  the  ASME  in  reference  25. 

The  flow  in  the  test  section  was  taken  to  be  equal  to  the  main  line 
flow  minus  the  flow  removed  by  boundary  layer  suction. 

The  free  stream  length  Reynolds  number  was  calculated  from  the  relation, 
Rex  = VxP,  where  |0  and were  taken  at  the  free  stream  temperature.  In 
the  casC  of  a heated  run  this  was  obtained  as  the  average  of  the  inlet  and 
exit  temperatures  of  the  test  section.  'Wien  calculating  the  mean  boundary 
layer  Reynolds  number  the  properties  were  selected  at  a temperature  equal 
to  the  average  between  the  wall  temperature  and  the  free  stream  temperature, 
the  wall  temperature  being  taken  as  the  average  of  the  wall  thermocouple 
readings.  For  the  friction  factor  data  the  mean  value  of  x between  two 
pressure  taps  in  the  laminar  zone  was  taken,  as  suggested  in  reference  6, 
as  the  average  of  the  arithmetic  and  geometric  means  of  the  axial  distances 
of  the  two  taps  from  the  boundary  layer  suction  slot  at  the  entrance  of  the 
test  section. 

The  local  apparent  friction  factor  between  two  pressure  taps  was 
obtained  from  the  formula, 

4f^pp  * Ap/A(x/D) 


where  Lp  is  the  static  pressure  change,  and  A(x/D)  is  the  number  of 
diameters  between  the  two  taps. 

The  values  of  PrpQTAL  — - were  obtained  directly  from  the 

manometers  as  the  difference  of  two  readings. 
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APPENDIX  B 


Approximate  Calculation  of  Experimental  Curve  Obtained 
by  the  Boundary  Layer  Probe  In  the  *Laminar  2one 

To  obtain  an  approximate  idea  of  the  shape  of  curve  obtained  by  the 
boundarv  layer  probe  in  the  laminar  zone,  a velocity  profile  for  flat  plate 
flow  was  selected  from  Schlichting^®^  (page  191).  Since  near  the  entry  of 
the  test  section  the  boundary  layer  is  thin,  a flat  plate  approximation, 
neglecting  the  effects  of  cross-section  curvature  and  pressure  gradients 
along  the  flow,  was  sufficiently  accurate  for  this  computation. 


The  velocity  profile  selected  was 


3IS 


In  the  case  of  flow  in  a pipe  entry,  XJ^  = U,  and  the  value  of  TJ 
changes  as  the  flow  develops  along  the  pipe  length.  In  the  present  example, 
since  the  boundary  layer  is  thin,  U can  be  considered  as  a first  approxima- 
tion equal  to  V,  the  mean  velocity  in  the  tube.  Then  (2)  reduces  to 


k(x)  - where  ^ex  - (3) 

Since  the  dynamic  pressure  at  a fixed  distance  from  the  wall  is 
proportional  to  u2,  we  obtain  from  (1), 


Substituting  (3)  into  (4)  results  in 
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Simplification  and  rearrangement  then  produces  the  final  form 


S1(<1oo  I.SSo.ooo 


Numerical  Example 

Since  the  experimental  probe  used  had  a finite  width,  the  dynamic 
pressure  indicated  was  actually  an  average  value  over  a certain  thickness 
of  the  boundary  layer.  To  simplify  the  calculations  a fixed  value  of 
y ~ 0.020  inches  was  used,  which  is  in  the  neighborhood  of  the  actual  probe 
size.  Other  typical  values  assumed  were: 

ReD  = 100,000  D = 2"  T0  = 75°P 

and 

V = — 5—  = 21,400  ft/hr 
or  V2  = 35.3  ft2/ sec2 

p 

Then  from  (6),  u was  calculated  at  each  value  of  x,  for  which  the  values 
PtotAL  — ^STATIC  a u2/2S  were  found  at  each  point  where  Rex  s Re^x/D). 

The  final  results  are  shown  in  Figure  6,  where  the  values  of 
PtotAL  ““  ^STATIC  ^ave  been  plotted  against  Re^. 
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APPENDIX  C 


Determination  of  an  Approximate  Expression  for  4?app  in  the 
Laminar  Entry  of  a Tube  When  the  Boundary  Layer  Thickness 
is  SmaTl~~Compared  With  the  Tube  Radius 


In  this  section  an  expression  for  4?^pp  In  the  laminar  entrance  of  a 

2 -3 

tube  is  derived  which  applies  for  values  of  = Re^/Re^  less  than  10  • A 

form  of  the  continuity  equation  is  first  obtained  by  considering  the 
physical  situation  which  exists  in  the  laminar  entry*  The  equation  is  then 
integrated  by  using  Pohlhausen’s  approximate  velocity  profile,  and  put  into 
a form  where  V/U  is  a function  of  SVd  and  dU/dx.  Successive  approximations 
of  U as  a function  of  x are  obtained  from  this  relation  and  are  used  with 
Thwaites'^^  approximate  method  of  laminar  boundary  layer  calculation  to 
obtain  higher  order  approximations.  The  calculations  proceed  as  follows: 
Continuity  Relation 

Tube  Wall 

Boundary 
Layer 


Assuming  steady  flow  and  applying  the  continuity  equation  for  a control 


volume  between  sections  1 and  2,  we  obtain  the  expression, 

x rl 

z 1>*V  « t{V-ZI)  \J  V TT  I (D-2^tcclj 


Simplifying,  we  get 


V.  - \ - 4/— \ t-  4/i\l  +-  f -_L_  f sw-As 
\j  ’ M>)  \d)  \)\)  \ UDM  j 

• o 


This  is  equivalent  to 


1.  I - 4 [ X Am  ♦ A.  f S is  t .!_  [ u-is-  _i_  f 
l)  L D 5 O'l.  UH\.  J VO').3  1 
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Regrouping  of  the  integrals  results  in  the  expression 

4 


V 

U 


Introducing  the  definition  of  displacement  thickness, 

>■  • (>  - 5) 


we  obtain 


1 = i - 4L*  v a. 


Me -?)’*> 


V 

To  evaluate  the  last  term  in  (2)  a velocity  distribution  u/U  as  a 


(2) 


function  of  y must  be  assumed.  A good  assumption  is  a polynomial  of  the 


fourth  order  due  to  Pohlhausen 


(18) 


_U 

u 


/i\  . ill)1 

UJ  UJ 


where 


Z"  dlx 


and 


+ 1(1)*  Jl\* 

UJ  UJ 


I « 1.  . 2L 

£ 10  110 


(3a) 
(3b, c) 


Substituting  (3a)  into  (2),  and  integrating  the  latter,  we  obtain 

V_  a \ - 4i*  + 1 £ . _X  £ 

U D 15  D1  45  0l 

Then,  from  (3c),  solving  for  l'  in  terms  of  and  X,  we  get 


(4) 


%%  - loo  p , too  S**  4.  so  » IS  + 

*1  SI  3(.m 


(5) 


If  the  expressions  (5)  and  (3b)  for  and  A are  now  inserted  into  the 
integrated  continuity  equation  (4),  the  latter  can  be  reduced  to  the  form 


V_  1-41*  + 2ooo(i!\Xl!liy  *... 

\J  D 21  VD  1 2U1  \ 0 Jr  1% 


(6) 


1st  Iteration 

To  begin  the  iteration  process,  an  expression  for  the  core  velocity 
U as  a function  of  x is  needed.  As  a first  approximation  the  potential 
flow  solution  for  flow  in  a pipe  is  assumed,  that  is  U = V,  or  V/U  = 1. 
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It  Is  noted  that  this  is  merely  the  first  order  term  in  the  continuity 
relation  (6),  If  we  substitute  this  into  Thwaites'  expression  for  the 


square  of  the  momentum  thickness 

el  * U5dU 

we  have 

0l  ^ ,45  yx 

V 

Continuing  with  Thwaites'  method. 


VM.  - - \J  0 - O 

Also,  from  Thwaites'  method 


as  U*  s o for  this 
1st  approximation 


%*  = 0 


(7) 

(3) 


(9) 


(10) 


By  plotting  H vs.  m as  given  in  table  I in  reference  21,  the  relationship 
H = H(m)  can  be  approximated  in  the  region  of  decreasing  pressure  by  the 
straight  line, 

H(m)  = 2.48m  +•  2.56  (11) 


Then  putting  (3), 


(9),  and  (11)  into  (10)  we  have, 

.4S*OC  ',l  ^.5U)  * I.U  JL. 

V I $7* 


(12) 


which  is  the  final  result  of  the  1st  iteration.  It  is  noted  that  this 
expression  is  the  same  as  that  obtained  in  the  Blasius  solution  for  flow 
over  a flat  plate.  This  might  have  been  anticipate^  as  very  close  to  the 
entry  of  the  pipe  the  boundary  layer  is  so  thin  compared  with  the  tube 
radius  that  it  approximates  flat  plate  behavior. 

2nd  Iteration 

We  now  take  as  a closer  approximation  the  first  two  terms  from  the 
continuity  relationship  (6), 


• or  U « i + «!V 

UP  V D 
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Substituting  %*  from  (12)  into  (6)  we  have 


U_  _ | f to.%%  Vd. 

V " ^ 


\ + to.n 


(13) 


This  expression  is  then  inserted  into  equation  (7),  thus  yielding 

A* 

Expanding  the  terms  in  square  brackets  by  means  of  the  binomial  theorem, 
and  integrating,  we  get 


& - .45 

1 + (..**f&T* 

-I*  r x 

i f (.,**{&! 

Ee*  . 

v l 

l £«B  . 

e"  - .45 

l - 41.1*  v.  1 

is 

r 

X + U.p  X +... 

V 

or 


S'  - .45  V*. f i - |15E  + ... 

v l e.„ 

From  (9),  (13),  and  (14)  we  have 

. [-5.44 

= - U 9 , l x g«p]L  V v fe0jJ 


(14) 


m. 


or 


y*. 


= -1.55  f | . |V55  + 


Ee« 


Ket 


From  (10)  and  (11), 


i*  * 0 (2.1?)^+  2.s0 


(15) 


(16) 


Substituting  (14)  and  (15)  into  (16)  and  retaining  only  terms  up  to  order 
fRe^/Rep,  there  is  obtained 


'b*  = 

. 4*  v* 

7i 

1 - ri*{i*< 

r 

2.4* 

_i.ss  f i_ 

+ 2.SL 

. V . 

to 

Ceo 

5 - 

or  . | 


nix 

(£7 


I - 10.U 


This  is  the  final  result  of  the  2nd  iteration. 


(17) 
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3rd  Iteration 


V/e  now  take  as  a closer  aporoximation  the  first  three  terras  frora  the 
continuity  relation  (6),  and  proceed  in  the  same  manner  as  before: 


V . I - 4V*  , m(1_1 

u T>  VI>J 

Substituting  (17)  for  Into  (19),  and  simplifying,  we  obtain 

V - I - (o.M )| Eex  f il.|  Eex  ^ 

U W 


(18) 


or 


(19) 


U _ \ » \o%%  l|  _ 4^1  ?e»  ^ . . . 

V ' ^ 

which  is  the  same  as  (13),  except  that  we  have  retained  the  next  higher 
order  term. 

Equation  19  is  now  Inserted  into  (7)  with  the  result,  after  integration, 
simplification,  and  dropping  of  all  terms  of  order  higher  than  Rex/Rep^,  that 


* {fZ 


I - Ifc.ss  )jFex_  + Eex 


Pe1 


(20) 


From  Eqs.  (9),  (19),.  and  (20),  we  have 


•-  - ,45 


S.44  1^.  _ ioU-1  Ket 


k,  it 

Solving  for  from  Eqs.  (10),  (11),  (20),  and  (21),  we  have,  after 
simplification, 

**  * 1-12  X 


(21) 


I - 10.  U ffoL  + M 

(feT  t ?*»  ^ 

Now  we  insert  (22)  and  dU/dx  from  (19)  into  the  first  four  terms  of  the 
continuity  equation  (6),  and  thus  we  obtain 

’ V.  , \ - 4i*  . ito/lt  t tooo  /|?jl  Ay 


(22) 


pi  Vd 


mi 


\J  * 

with  the  result  that 

v . u.n(tT.  » io.i  . mi. 


U 


f=e.' 


fe,1 
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or 


VL  t l + . 43.  S Rex  4 loUo  Rex^1  _ . ..  (oti 

v w V 

It  Is  noticed  that  the  successive  coefficients  in  the  series  grow 
quite  rapidly  and  thus  good  convergence  is  obtained  only  in  the  region  where 
small  and  less  than  10“^.  A similar  expression  was  obtained 
by  Atkinson  and  Goldstein  (see  reference  5,  pages  305-306),  their  result 
being 


. »t. sis  4 iiss.SU  Rt«v' 

V " ft* 

We  may  now  write  the  Bernoulli  equation  for  the  frictionless  core  of 
the  flow  between  sections  1 and  2 in  the  figure  on  page  36.  Thus, 


K-  h 
i-?v' 


where  and  pg  are  the  static  pressures  across  the  whole  cross-section 
at  1 and  2 respectively.  Therefore,  from  the  definition  of  4?App, 


4 A.  * |b\‘  - I 

* t>  \\J 


(24) 


Putting  (23)  into  (24),  we  have  the  final  result. 


l - 2.Vi  + ui  ... 

(25) 


A comparison  of  points  calculated  from  this  relation  with  the 
experimental  curve  obtained  by  Kline  and  Shapiro^® ) is  shown  in  Figure  30. 


I 
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APPENDIX  D 


Friction  Factor  in  the  Laminar  Inlet 
of  a Smooth "Tube  Up  to  the  Point  of  Formation 
o7  Polseullle  Flow 

Form  of  the  Boundary  Layer  Velocity  Profile 

As  discussed  previously  in  section  VIII,  the  determination  of  the  growth 
of  a laminar  boundary  layer  in  the  entry  of  a tube  is  complicated  by  the  fact 
that  the  boundary  layer  growth  affects  the  velocity  in  the  frictionless  core 
of  the  flow.  The  growth  of  the  layer  produces  an  acceleration  of  the  core 
velocity,  thus  causing  a pressure  gradient  to  exist  in  the  flow  direction. 
This  in  turn  affects  the  boundary  layer  growth.  Since,  for  a given  flow, 
the  core  velocity  (and  hence  the  pressure)  at  a particular  position  along 
the  tube  is  dependent  on  the  thickness  of  the  boundary  layer  in  relation 
to  the  pipe  radius,  we  might  expect  that  the  boundary  layer  velocity  profile 
would  also  depend  upon  this  ratio.  Thus  in  the  case  of  pipe  flow  we  might 
say  that. 


= Function 


(Ml 


We  may  then  assume  that  u/U  is  of  the  form, 

4 d.|  j + • • * 4 +■  +■  * • • 

where  aQ,  a^,  b-^,  bg,  c^,  c-^g,  etc,,  are  constants.  For  the  sake  of 
simplicity  we  shall  retain  only  terms  up  through  the  fourth  order,  and 
therefore  we  choose  a polynomial  for  the  velocity  profile  in  the  form 


(1) 


Ji  - a„  + a, 
l) 


S)  “'®'  • “S'  • MS'  * 48  • Mff  * MS  • MW 


m '-©18  • MS®  • '■ffi'B 

where  the  constants  have  to  be  fixed  to  satisfy  the  boundary  conditions 


(2) 


of  the  problem. 
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Boundary  Conditions 


The  first  boundary  condition  to  be  considered  is  that,  at  the  wall, 
the  fluid  velocity  must  equal  zero.  Thus,  aty=0:  u = 0,  Inserting 
this  condition  into  (2)  we  have 


A second  boundary  condition  which  we  shall  impose  is  that  the  velocity 
u at  the  edge  of  the  boundary  layer  shall  be  equal  to  the  core  velocity  U. 
Thus,  at  y s £ : u * U,  Inserting  this  into  (4)  we  have 

\ a 0-,  + + a.*  v (l\  v C.  ,xfL\X  v v txi  (L\  f + C wi  /S\  (5) 

\%)  UJ  \%)  \v  W Vr) 


and  we  obtain  the  following  equations  by  equating  terms  containing  equal 
powers  of  (&/R): 


a,  + a,  + 

+ =1 

(6) 

C„  + Ci,  f 

* 0 

(7) 

Cu  v Cix 

= O 

(3) 

*■  O 

(9) 

The  third  boundary  condition  imposed  is  that,  at  y *k:  ^u/iy  s 0. 

Dropping  the  C13  term  according  to  (9),  and  differentiating  (4),  we  have 


Applying  the  boundary  condition  leads  to 


0 * 


a,  f i(Kx 


♦ 3a5  + 4a.<  v Cu 


ft  + 


(11) 
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which  provides  the  additional  equations 


a,  + lOii  + 3cl5 

* o 

(12) 

C„  + 2Ct,  + 5Csi 

• o 

(13) 

ClX  + ^Cxx 

• o 

(14) 

From  (3)  and  (14)  we  find  that 

c12  an<^  c22  mus^  ©ach  equal  zero. 

Following  the  entry  region,  within  which  the  boundary  layer  thickness 
increases  and  eventually  fills  the  tube,  the  velocity  profile  acquires  the 
parabolic  shape  characteristic  of  fully-developed  laminar  flow.  Thus  we 
shall  impose  a fourth  boundary  condition,  which  states  that  when  ^ = R, 
the  velocity  profile  assumes  the  parabolic  shape  given  by 


(15) 


This  condition  is  then  applied  to  equation  (4)  (from  which  the  terms 
involving  0^3,  C12*  an<*  c22»  have  already  been  omitted  as  they  have  been 
found  equal  to  zero)  and  we  obtain 

* . m/I\*  * . t„  h\ 


U)  UJ  u 

+ Cxi  [j}.\  V 


w 


b 


u 


Thus  we  have  the  equations 


(1\S 

UJ 


(16) 


OL,  + C„  =2  (17) 

Q-t  + Cxi  - - 1 ( 13 ) 

fl-i  +CV  ■ O (19) 

*0  (20) 

Since  a4  is  equal  to  zero  we  are  finally  left  with  a profile  of  the  fonn, 

- a,|I\  v o-i[.L\  + \ t (21) 

\j  UJ  UJ  UJ  \&JUJ  UJ  UJ  UJ  w 

which  has  six  unknown  constants.  It  would  seem  that  these  could  be  solved 


for  from  the  seven  equations  (6),  (7),  (12),  (13),  (17),  (18),  and  (19). 


These  become,  after  dropping  the  constants  which  are  equal  to  zero, 

OL,  + * A-*  * I 


(22) 
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Ci,  4 Cxi  4 Csi  a O ( 23) 

Q-i  4 Zdx  4 * 0 ( 24 ) 

C„  4 £Cn  4 3C^|  - o (25) 

0.,  4 Ci,  -1  (26) 

C lx  +C»,  *--|  (27) 

4?j|  *0  (23) 


However,  if  (26),  (27),  and  (23),  are  inserted  into  (22),  equation  (23)  is 
obtained.  Also  if  (26),  (27),  and  (28),  are  inserted  into  (24),  then  (25) 
is  obtained,  so  that  only  the  five  equations  (22),  (24),  (26),  (27),  and 
(28)  are  independent.  Thus  one  more  equation  is  needed  before  the  constants 
can  be  determined.  This  will  be  found  from  a fifth  boundary  condition 
derived  from  the  steady-state  equation  of  motion  in  cylindrical  coordinates. 


br  V**1  r ***  bxv 

The  boundary  condition  states  that,  at  the  wall,  where  y = 0:  u * 0 and 


v a 0.  We  neglect  the  term  in  the  equation  as  is  usual  in  boundary 
layer  theory,  and  apply  the  Bernoulli  equation  in  the  frictionless  core  of 


the  flo?/  to  evaluate  bp/bx.  Thus  we  have. 


0 * It  + 4 

^ bx  V "Vr1  r V / 

with  - L tt  a ^ All  and  fc.  *■  C 4 M 

fi  ^x  lx.  J 

Substituting  into  (30)  and  keeping  in  mind  that 
y s Owe  must  fulfill  the  condition. 


(30) 

where  R is  the  pipe  radius. 

^ x - , we  have  that  at 

■by 


tfu  4 _L_  It  x - t it 

We  now  proceed  to  evaluate  the  terms  in  this  equation, 
the  velocity  profile  (21)  we  have, 


(31) 

Differentiating 


_ y_ 
^ ' b 


OL,  4 20.x(l\  4 51 ifctf  4 c„/i\  4 'Cv,(I\(V\  4 iCvhV/d 


i\W 


lb)  Ip; 


(32) 


m u 

cT 

L 

Ui. 

T°  * 

[ WJ 

(33) 
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and 


. u 

t + (»&*(— \ + 2C 

'T 

s>i 

1 

\\J 

u 

r 

lAt  V 1Ci\fl) 

W\ 

3ao 

i WJ 

(34) 

(35) 


It  is  next  necessary  to  evaluate  the  terra  dU/dx,  and  this  can  be  obtained 
from  the  continuity  relationship  (Eq,  (1),  Appendix  C).  This  can  be 
written  in  the  form, 


V _ \ - i|s]\ 
U U) 


r*  r ^ 

11  - 2_  I it  h 

*1 y p*L » 3 


(36) 


After  evaluation  of  the  integrals  the  equation  is  solved  for  U,  the  result 
be ing  / 

U*  V/  + (3?) 

- i - % - 1 x'){$ + (-  t «•  - * 1 c*)  ftf] 

This  can  then  be  differentiated,  and  dU/dx  is  found  in  the  form. 


AL  (33) 

dx 

This  then  presents  the  problem  that  di/dx  must  be  evaluated,  but  once  this 
is  found  the  entire  problem  is  solved  after  the  constants  a^,  c-^,  etc., 
have  been  evaluated.  To  find  dS/dx  we  must  turn  to  the  momentum  equation, 
which  may  be  derived  for  a tube  in  momentum- integral  form  by  considering 


x 

d.x 


a control  volume  as  shown  in  the  figure  below. 

Control  Volume — ^ . ^Boundary  Layer 


Tube  Woll- 


If  the  net  force  on  the  control  volume  is  equated  to  the  net  outgoing 
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flux  of  momentum,  we  obtain  the  equation 


- ma - * dil  puv<uu  di 

^ icy  io 


(39) 


Applying  the  Bernoulli  equation  in  the  frictionless  core  to  evaluate 


dp/dx,  and  using  X 


U <JIU  . 2v/i  U.' 


, (39)  can  be  reduced  to  the  form 


<U  Z v»y) 


“ Tip), 


(40) 


All  of  the  terms  in  this  equation  can  then  be  evaluated  by  the  use  of  the 
assumed  velocity  profile  (21),  The  first  term  on  the  left  hand  side  and 
both  terms  on  the  right  will  involve  the  factor  d&/<ix  since  these  terms 
contain  differentiations  with  respect  to  x.  The  equation  can  then  be 
solved  for  db/dx  and  its  value  inserted  into  (38),  Equations  (33),  (35), 
(37),  and  (38)  are  then  substituted  into  the  boundary  condition  (31)  and  y 
is  set  equal  to  zero.  The  condition  is  then  examined  to  attempt  to 
determine  the  relationship  between  the  undetermined  constants  so  that  they 
may  be  evaluated.  This  procedure  is  extremely  long  and  involves  very  large 
algebraic  expressions,  so  it  will  be  omitted  from  this  discussion.  The 
results  will  be  stated  and  these  results  will  be  later  verified  by 
substitution  into  the  boundary  condition. 

It  was  found  impossible  to  fully  satisfy  this  boundary  condition,  so 
it  was  decided  to  satisfy  it  partially.  The  calculations  showed  that*  the 
left  hand  side  of  (31)  had  a constant  term  a2,  with  all  other  terms 
multiplied  by  (S/R)  or  (S/R)^,  On  the  right  side,  however,  no  constant 
term  appeared  and  all  terms  were  multiplied  by  (S/R)  to  the  first  or  a 
higher  Dower,  Thus  in  order  to  partially  satisfy  the  boundary  condition, 
ag  was  taken  equal  to  zero.  This  means  that  the  boundary  condition  was 
completely  satisfied  in  the  region  where  (S/R)  is  sufficiently  small  that 
all  terms  except  the  constant  term  may  be  neglected.  The  boundary  condition 
was  partially  satisfied  elsewhere. 

With  a2  equal  to  zero,  the  constants  a^,  a^,  c^,  Cgp  and  c^  can  be 
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evaluated  from  equations  (22),  (24),  (26),  (27),  and  (28),  resulting  in 
the  values 


a^  s 3/2,  * 8,  a<j  * **1/2,  ^11  ™ 1/2,  ^21  " *1 , and  ™ 1/2  • 


Substituting  these  into  equation  (21),  we  have  as  a final  result 

\)  7 v%)  ns 


zUAPi  VbMPr  zVSjVP 


(41) 


To  interpret  this  expression,  the  terms  may  be  rearranged  in  the 


form 


(42) 


Prom  this  it  is  seen  that  as  (S/R)  goes  to  zero,  u/U  becomes  equal  to 


IM1- 


This  is  a familar  cubic  velocity  profile  for  a flat  plate 
(see  reference  18,  page  191),  and  thus  the  profile  approaches  flat  plate 
behavior  when  the  boundary  layer  is  very  thin  in  comparison  with  the  pipe 
radius.  As  (S/R)  goes  to  one,  the  first  term  on  the  right  side  goes  to  zero 
while  the  second  term  increases  in  magnitude.  When  (S/R)  is  equal  to  one 
the  profile  becomes  u/U  * which  is  the  parabolic  profile  for 

fully- developed  laminar  flow.  Thus  as  ($>/R)  g°es  from  zero  to  one»  the 
shape  of  the  profile  gradually  changes  from  a cubic  flat  plate  profile  to 
a fully-developed  parabolic  profile. 


Evaluation  of  Terms  in  Continuity  and  Momentum  Equations 

Using  the  approximate  velocity  profile  (41)  we  will  now  proceed  to 
solve  the  continuity  and  momentum  equations  to  determine  the  growth  of  the 
boundary  layer  as  a function  of  distance  along  the  tube.  Substituting  the 
values  of  the  constants  into  (37),  which  has  been  obtained  from  the 
continuity  relation,  we  have 


u - 


Differentiation  of  this  expression  results  in 

AU  . .V  u 


» 


(44) 
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No?/,  for  convenience  let 


(45) 

(46) 


The  first  terra  in  the  momentum  equation  (40)  can  then  be  expressed  as, 

U AO  . - v.1  &.  ii  (47) 

A*  K-  6'  A* 

To  evaluate  the  second  term  in  the  momentum  equation,  we  determine 
f^\  from  equation  (33),  with  the  result  that 


Now  let 


iu\  - fl  V 1 (±\(  V 

Up)] 

"(1*4©] 

so  that  the  second  term  in  the  momentum  equation  (40)  may  be  written 

- - - \y  V _E 


(48) 

(49) 


£ 


n t 


(50) 


Now  consider  the  third  term  in  the  momentum  equation. 


z_  L 


VfV  hHj 


(51) 


This  is  differentiated  by  parts  to  obtain 

* T C rl . _ \ 


\)'L 

A% 


h ’ i (sin*’ 


'*  I L 

Evaluating  the  integral,  we  have 

l 

kl  /V\  - loo^  / 

2.5  VpJ  ILloi 


z\j  in 


(52) 


if  - 11 1!]  _J_/i 
kpJ  W«VW  SUoVPi 


(53) 


tW^'' 

Substituting  this  into  (52),  introduping  U anc$  dU/dx  from  continuity,  and 
I ‘ 

carrying  out  the  differentiation, we  obtain 


* i 1 C©>  •»  *1]  ■ [Ui  { a - *§!  © - & - k (8TV 


as 

77 


SjfU  _ 1003  /s\l.  il  (if  _ J_  (if 
is  VP/  KAoVW  %1oVP/  SUoW 


-li'i 


AS 

(2.  S1]  A* 


(54) 


Now,  for  convenience,  let 


(55) 

(56) 


so  that  the  third  term  in  the  momentum  equation  becomes 

a 


£ 

r ■ 

c 

n 

IV1  k ’ 

U\ 

.*] 

Ax 

L R V. 

or 


ii 

Ax 


£) 

r 


C.  . ZGA 
6l  B* 


(57) 


The  fourth  term  of  (40)  is  then  considered, 

A 

in 

Carrying  out  the  differentiation,  we  have 


f.-ittUttf] 

o1 

w W 1 

v w 

(58) 


l- 


llQ  *(1\X 

[£)  VP-) 


z\)4u  vU  . i.ifi: 
dLx  R I VR. 


ii 

Ax 


(59) 


Substituting  into  this  the  values  for  U and  dU/dx  from  continuity  results  in 

(60) 


ii 

Ax 


‘-2^  A ’ 

a 

fl  1 

\W  \pj  J 

R S\ 

Ax 

R 

. V. 

WJ 

Now,  setting 


F» 


-l  + l[%} 

VR> 


(61), (62) 


this  term  becomes 


a 

[t 

P-iA©  * F " 

c. . 

(63) 


By  equating  the  four  terms  just  evaluated,  the  momentum  equation  can 
be  written  in  the  form,  ^ 


Rearranging,  and  letting  Re^  ~ VD/V» 

-i_  ^ ii  [A-  + - - 1£A  - + 1 


?er 


U & 


v-  e> 


A 

It 


(65) 


Now  let 


3?li\  s fi_  + L - 1M  - 1A& 

^vpJ'U1  & bv  bx  bJUj 


t 


(66) 


so  that  we  obtain 


ft  n 
w 


_ A_  . 

Separating  variables,  we  may  now  integrate  the  expression 

■W*  . ,N 

fft  4ft  - 


(67) 


(68) 


and  find  Rex/ReD2  as  a function  of  (S/R).  This  has  been  done  graphically 
by  calculating  at  several  values  of  (S/R).  A curve  of  & vs«  (S/R) 
is  then  plotted,  and  the  area  under  the  curve  up  to  any  (&/R)  is  the  value 
of  -8Rex/ReD2  at  that  (S/R).  Also,  from  continuity  (43),  we  may  find  the 
value  of  U/V  at  that  (S/R)*  Then  from  Eq.  (24),  Appendix  C,  we  have 


(69) 


Thus  a value  of  4?^pp(x/D)  is  obtained  at  each  value  of  Rex/Rep  , and  a 
curve  may  be  plotted  which  is  shown  in  Figure  31.  Numerical  values 
obtained  during  the  graphical  evaluation  are  tabulated  in  Table  III. 
Verification  of  the  Fifth  Boundary  Condition 

We  shall  now  return  to  the  last  boundary  condition  and  cive  a 
verification  that  the  constant  ag  is  equal  to  zero  as  previously  stated. 
The  boundary  condition  to  be  fulfilled  was  that  at  y = 0: 


£±  +j_  h * -U. 

j-R  vj  ^ d* 


(31) 


Substituting 
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TABLE  III 


Summary  of  Factors  Obtained  In  Graphical  Integration 


(S/R) 

A 

B 

C 

£ E 

F 

G 

0 

-.75000 

1.00000 

0.97143  1. 

00000  1.50000 

-2.00000 

0 

0.005 

-.74717 

0.99625 

0.96546  0. 

99003  1.50250 

-1.99000 

.00484 

0.010 

-.74434 

0.99253 

0.95948  0. 

98010  1.50500 

-1.98000 

.00965 

0.050 

-.72192 

0,96320 

0.91143  0. 

90250  1.52500 

-1.90000 

.04707 

0.100 

-.69433 

0,9  2 780 

0.85084  0. 

81000  1.55000 

-1.80000 

.09114 

0.300 

-.58900 

0.79960 

0.60242  0. 

49000  1.65000 

-1.40000 

.23662 

0.500 

-.49167 

0.69163 

0.34405  0.25000  1.75000 

-1.00000 

.33143 

0.600 

-.44600 

0.64480 

0.21103  0. 

16000  1.80000 

-0.90000 

.35921 

0.800 

-.36067 

0.56427 

-.06290  0. 

04000  1.90000 

-0.40000 

.37420 

1.000 

-.28333 

0.50000 

-.34762 

0 2.00000 

0 

.33333 

(VR) 

* 

(&/R) 

Re^/Re^  4?^pp(x/D) 

0 

0 

0.050 

0.0000318 

0.07736 

0.005 

-.00094 

0.100 

0.000135 

0.16169 

0.010 

-.00192 

0.200 

0.000638 

0.34874 

0.050 

-.01071 

0.300 

0.00164 

0.56406 

0.100 

-.02425 

0.400 

0.00328 

0.81046 

0.300 

-.10544 

0.500 

0.00566 

1.09049 

0.500 

-.22315 

0.600 

0.00887 

1.40519 

0.600 

-.29073 

0.700 

0.01299 

1.75569 

0.800 

-.42721 

0.800 

0.01794 

2.14073 

1.000 

-.53651 

0.900 

0.02364 

2.55730 

1.000 

0.03002 

3.00000 

v;e  have. 
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Substituting 


we  obtain. 


jy  a _ v 

<£.x  ?. 


1l 

B' 


li 

ckx 


_U  A(J 

^ Tx 


But, 


where  Re^  = V%/p, 


Then  (73) 


reduces  to. 


(72) 


(73) 


T_  | X\  v - (— \ (74) 

4 U)  4W  W ^ vcft  “tAfc  vFfc] 

If  a2  had  not  been  equal  to  zero  there  would  be  a constant  term  on  the  left 
hand  side  of  (74).  However,  there  are  no  constant  terms  on  the  right  hand 
side  as  the  first  term  will  be  of  order  (&/R).  Thus  to  partially  satisfy 
the  relation  (74),  ag  had  to  equal  zero  as  stated  previously. 

Continuation  of  Curve  of  4T^pp(x/D)  vs,  Re^/Rep2 

On  Figure  31  the  boundary  layers  meet  when  4?^pp(x/D)  is  exactly  equal 
to  three,  at  which  point  (x/D)  » 0.03002  Rep  (the  latter  value  is  a little 
larger  than  the  value,  (x/D)  ■ 0.02875  Rep,  which  was  obtained  by 
Schiller^® ) ) . Beyond  the  point  where  (£/R)  = 1 we  may  extend  the  curve 
by  using  the  Hagen-Poiseuille  Law.  Thus  we  write, 

Mffcpp  A * At.  - -•  !*i  Lks  . U4bf&x\  (75) 

D i^v'  Pwwwufc  D feo  ?{,  \\fcpx) 
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where  A signifies  that  this  is  the  additional  friction  factor  to  be  added 
to  47^pp(x/D)  s 3.  Thus  for  example  when  * 0.07;  &47App(x/D)  = 

64(0.07)  = 4.48.  Then  47App(x/D)  = 3 + A47App(x/D)  = 3 4 4.48  = 7.48. 

Thus  a ooint  is  plotted  at  Re^/Re-^  - 0.07  4 0.03  a 0.10  and  47App(x/D)  = 
7.48,  through  which  the  curve  may  be  extended.  The  theoretical  curve 
obtained  (Figure  31)  follows  the  data  of  Kline  and  Shapiro^  for  the 
beginning  of  the  entrv  region,  and  then  gradually  approaches  the  Hagen- 
Poiseuille  line. 

Behavior  of  the  Solution  as  (S /R)  Approaches  Zero 

To  determine  how  the  solution  behaves  in  the  region  very  close  to 
the  entry  of  the  tube  as  (&/R)  goes  to  zero,  we  begin  with  expression  (66) 
and  neglect  terms  of  order  higher  than  (S/R).  This  results  in. 


Integrating  according  to  (68)  we  have, 

rVK 

-2L  (i\  ilk]  s -t 

1 \r)  W Zt* 

* o 

which  becomes, 


(76) 


?«.  S J5_  (1' 

> 1120  VR 


1 


(77) 
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Thus, 


(!\  s 1.21 


as  lk\ 


0 


(73) 


From  continuity  (43)  we  obtain, 

u . 


V 


Then  we  have, 


APP 


l 

- 1 1.21  ^ 

Q 

['-!©  ♦-]" 

\+  t-n.  ftZ 

- * Pf  - 1 

D Iv) 

* 13.12  JPe* 
^0 

(79) 


(BO) 


which  is  valid  near  the  pipe  entry.  This  is  shown  as  a dashed  line  on 
Figure  31,  and  may  be  compared  with  the  experimental  value  13. 74 fRe^/Re-p 


of  Kline  and  Shapiro 


(6) 


As  a check  on  this  result,  we  have  for  the  flat  plate  profile 


u/U  = U -p- 

*UJ  i\i 

w = 


(see  reference  13,  pages  190-191), 


£25. 

U„ 


2 *T  I 
31  V U. 
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\ 31_ 

I Z%0 


l\0  TSX 


13  \ M 


Vi 


(91) 


1.Z1  ffe. 

'W/y']  fe 

which  is  equivalent  to  (78), 

Local  Apparent  Friction  Factor  Along  the  Tube 

The  method  developed  has  led  directly  to  an  expression  for  the  mean 
apparent  friction  factor  7^pp,  but  it  is  also  of  interest  to  observe  the 
variation  in  the  local  apparent  friction  factor,  f^pp#  along  the  tube. 

This  can  be  found  by  differentiating  the  log-log  plot  for  mean  apparent 
friction  factor  in  the  following  manner.  Beginning  with  the  definition  for 
mean  aooarent  friction  factor  we  have, 

i rv» 


^kpp  * 


(*/  o) 

y 


^tapp  l- 


(82) 
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Differentiating  results  in. 


which  is  equivalent  to. 


By  measuring  the  slope  of  the  mean  apparent  friction  factor  curve 
(Figure  31),  it  is  oossible  to  determine  values  of  the  parameter  4f^ppRep 
at  various  values  of  Rex/ReD^.  The  universal  curve  thus  obtained  is  shown 
in  Figure  32,  and  is  in  agreement  with  the  experimental  data  of  Kline  and 
Shapiro (the  experimental  data  extend  up  to  Re^Re^  a 10~^). 

These  results  can  also  be  plotted  as  4f^pp  vs.  Rex.  In  this  case  a 
family  of  curves  is  obtained  (Figure  33),  each  curve  being  for  a different 
diameter  Reynolds  number. 

growth  of  the  Laminar  Boundary  Laver  '.Vi  thin  a Tube  Bn  try 

The  values  of  (S/R)  vs.  Re^Re-p^  which  were  obtained  from  the  graphical 
integration  process  have  been  plotted  in  Figure  34.  The  curve  illustrates 
the  manner  in  which  the  boundary  layer  thickness  increases  along  the  length 
of  the  tube.  For  comparison  purposes  the  flat  plate  relation  (81)  has  also 
been  plotted  on  the  same  figure.  The  curves  illustrate  that  the  favorable 
pressure  gradient,  which  is  present  as  the  flow  develops  within  a tube  entry, 
causes  the  boundary  layer  to  grow  at  a decreased  rate  as  compared  with  the 
growth  on  a flat  plate  with  no  external  pressure  gradient. 

Substitution  of  the  Pohlhausen  Profile  for  the  Cubic  Profile 

It  is  recalled  that  as  (&/R)  approaches  zero,  the  expression  (42) 

reduces  to  a cxibic  flat-plate  velocity  profile.  It  was  thought  that  a 
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better  approximation  might  be  obtained  if  the  fourth  order  Pohlhausen 
profile  were  substituted  for  the  cubic  expression.  Thus  a new  profile  was 
formulated, 


u 

u 


\i)  \b 


\-(» 

VRJ 


■fl-Bf 


(94) 


This  expression  was  used  to  evaluate  the  terms  in  the  continuity  and 
momentum  equations  in  a manner  similar  to  the  computations  previously 
carried  out.  Substitution  into  continuity  yielded  the  expression, 

U:  ! 

? 


where. 


?* 


\ - —( + _L  ( A\  v _L  ( ^ 
sW)  IS\rJ  SoIrJ 


(85) 

(86) 


Differentiating  with  respect  to  x gave, 

W ..  J Q U 
<U  t P'  A* 


where. 


Q * 


.i  + j_/i\  + ±IV\ 

S ISVr)  I<Ar 


(87) 


(88) 


The  first  term  in  the  momentum  equation  then  became, 

UiU.-V'Q.ii  (39) 

h E P’  h 

The  second  term  in  the  momentum  equation  could  be  immediately  evaluated  as 


-sylhA  . -ih  V (90) 

R uV'j.o  ft  P 

After  evaluation  of  the  integral,  the  third  term  in  the  momentum  equation 
assumed  the  form. 


where 


A 

V <u 


n 

Y 

T 

?sa 

dLx 

1*1 

?'  . 

(91) 
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(92) 
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Finally,  the  fourth  terra  in  the  momentum  equation  became. 


(93) 
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(Ax 

Equating  the  terms  and  simplifying  yields  the  momentum  equation  in  the 
form. 


(94) 


* ±£ 


Setting 


Ax 

K*t> 


0,4.1.  ?SQ  . 2Q.S  4-1 
pl  P ?'  P'  ? 


I 


(i\  - 20^  4 £ 


P' 


Pl 


P 


produces  the  equation, 

•£z(« 


*e0  Ax 

This  expression  can  then  be  integrated, 

;Vn 

-IU 


(95) 


(96) 


(97) 


(98) 


This  has  been  carried  out  graphically  as  in  the  previous  computation,  and 
the  principal  results  are  shown  in  Table  IV.  The  results  were  found  to 
lie  on  the  same  curve  (Figure  31)  which  was  obtained  using  expression  (41), 
except  that  the  entrance  length  was  slightly  smaller,  the  value  being 
x/D  s 0.02961  Re^,  as  compared  with  the  earlier  value,  0.03002  Re^. 

To  determine  the  behavior  of  this  solution  as  (S/R)  approaches  zero, 
we  follow  the  same  procedure  that  was  used  to  evaluate  expression  (78) 
with  the  result  in  this  case, 

dL  - n' 

Etc 

(86)  we  obta  , 


^ £ ll.W 


as 


(99) 


From  continuity  (85)  and  (86)  we  obtaifi, 
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TABLE  IV 


Summary  of  Factors  Obtained  in  Graphical  Integration 


(*/R) 

P 

Q 

S 

T 

Z 

0 

1.00000 

-.60000 

0 

1.16508 

0 

0.005 

0.99700 

-.59933 

0.00530 

1.15619 

-0.00119 

0.010 

0.99401 

-.59866 

0.01156 

1.14732 

-0.00242 

0.050 

0.97017 

-.59303 

0.05604 

1.07663 

-0.01354 

0.100 

0.94070 

-.58567 

0.10768 

0.98906 

-0.03092 

0.300 

0.82690 

-.55100 

0.27110 

0.64743 

-0.14679 

0.500 

0 . 72083 

-.50833 

0.36756 

0.31944 

-0.35705 

0.600 

0.67120 

-.43400 

0.39153 

0.16047 

-0.50526 

0.800 

0.57973 

-.42933 

0.39259 

-.14763 

-0.89555 

1.000 

0.50000 

-.36667 

0.33333 

-.44286 

-1.37460 

(*>/R) 

Rex/ Re  D2 

4T^pp(x/D) 

0.050 

0.000020 

0.06244 

0.100 

0.000089 

0.13005 

0.200 

0,000422 

0.28276 

0.300 

0.00111 

0.46249 

0.400 

0.00230 

0.67442 

0.500 

0.00414 

0.92455 

0.600 

0.00682 

1.21971 

0.700 

0.01052 

1,56739 

0.800 

0.01533 

1.97539 

0.900 

0.02172 

2.45105 

1.000 

0.02961 

3.00000 
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Then  we  have, 

i . |y\  . | ■£  H.oifre,  (loo) 

D W Pe. 

as  compared  with  the  value  13.92'|Rex/Rej)  obtained  in  the  earlier  calcula- 
tions (Eq.  (80)).  Equation  (99)  may  be  checked  by  the  use  of  the  constants 
listed  in  reference  18  (pages  190-191),  for  the  flat  plate  profile, 
u/U  a + • Thus  we  have. 


U,  i u. 


iii  rEx 
SL  U. 

3IS 


£\  - [HE*  SI  . 

*)  1 ^ 1T 


(101) 


The  boundary  layer  growth  for  these  computations  is  plotted  as  a 
set  of  dashed  lines  on  Figure  34* 

Correction  Factor  for  a Viscometer 

In  a viscometer  the  rate  of  flow  and  pressure  drop  through  a tube 
leading  from  a tank  is  measured,  and  from  this  the  fluid  viscosity  can 
be  determined.  The  calculations  are  based  on  the  assumption  that  fully- 
developed  Poiseuille  flow  exists  throughout  the  tube.  However,  due  to  the 
fact  that  the  flow  is  developing  in  the  entry  region  of  the  tube,  a 
correction  factor  must  be  applied,  as  the  actual  loss  in  this  region  will 
be  larger  than  that  calculated  from  the  fully- developed  laminar  equations. 
The  pressure  drop  equation  used  for  these  calculations  is  of  the  form. 


— - ^4  111  + l + Correction  Constant  (102) 

The  first  term  on  the  right  represents  the  pressure  loss  in  the  tube  if 
fully-developed  laminar  flow  were  to  exist  throughout  the  entire  tube 
length.  The  second  term  represents  the  change  in  kinetic  energy  from  the 
tank  to  the  tube  entry,  and  the  third  terra  is  the  correction  factor  which 
must  be  applied.  This  correction  factor  may  be  easily  calculated  by 
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taking  the  actual  total  friction  drop  in  the  entry  region  of  the  tube 
(up  to  (E/R)  s l)  and  subtracting  the  friction  drop  that  would  occur  in 
that  length  for  fully-developed  laminar  flow.  Thus  we  have  for  the  case 
when  velocity  profile  (41)  was  used, 

Correction  * 3.000  — 64(0.03002)  ■ 1.079 
For  the  velocity  profile  (84)  we  have. 

Correction  • 3.000  — 64(0.02961)  * 1.105 
Values  for  this  correction  have  been  calculated  by  other  investigators, 
their  results  being  as  follows:  Schiller ^ 16 ),  1.16;  Langhaar(7),  1.28; 

Boussinesq^1),  1.24;  and  Atkinson  and  Goldstein^ 5),  1.41.  Experimental 
values  summarized  by  Langhaar^7),  range  from  1.24  to  1.32,  so  it  apoears 
that  the  values  computed  by  the  present  method  may  be  slightly  low,  with 
the  Pohlhausen  type  approximation  (84)  producing  a slightly  better  result 
than  the  profile  of  equation  (41).  In  the  actual  experiments  there  may 
be  a frictional  loss  at  the  tube  entry,  which  might  account  partially 
for  the  difference  between  the  observed  values  and  the  present 
computations. 


FIG.  I SCHEMATIC  DIAGRAM  OF  THE  FLOW  SYSTEM 


FIG.  2 THE  STILLING  CHAMBER.  THE  FLOW  WAS  STRAIGHTENED  BY  A 
HONEYCOMB  CONSISTING  OF  TIGHTLY  PACKED  EXTRUDED  ALUMINUM 
TUBES.  THIS  WAS  FOLLOWED  BY  15  STAINLESS  STEEL  TURBULENCE 
DAMPING  SCREENS  (2  SHOWN  IN  FIGURE)  WHICH  WERE  EACH  SOLDERED 

TO  AM  IMTPPl  nririMA  AM  fell  II  AD  ObAuve  bima 


FIG.  3 DOUBLE  CONTRACTION  NOZZLES  AND  BOUNDARY  LAYER  SUCTION 
AT  THE  TEST  SECTION  ENTRY. 


FIG.  4 TEST  SECTION  ENTRY  SHOWING  THE  CONFIGURATION  OF  THE 
BOUNDARY  LAYER  SUCTION  SLOT  AND  PRESSURE  TAPS. 


FIG.  5 BOUNDARY  LAYER  PROBE.  THIS  WAS  INSERTED  THROUGH  THE 
EXIT  END  OF  THE  TEST  SECTION.  IT  MEASURED  THE  TOTAL 
PRESSURE  IN  THE  BOUNDARY  LAYER  AT  A FIXED  DISTANCE 
FROM  THE  WALL,  AS  WELL  AS  THE  STATIC  PRESSURE  IN  THE 
CORE  OF  THE  FLOW. 


FIG.  6 CURVE  CALCULATED  FOR  LAMINAR  ZONE  AS  INDICATED  IN 
APPENDIX  B.  TOTAL  PRESSURE  PROBE  ASSUMED  FIXED 
0020“  FROM  WALL.  Ri0 « 100,000  T0  • 75#F 


LOCAL  APPARENT  FRICTION  FACTOR,  4f 


FIG. 7 ADIABATIC  RUN  WITH  TRANSITION  DETERMINED  FROM  LOCAL 
APPARENT  FRICTION  FACTOR  MEASUREMENTS.  TO  BE 
COMPARED  WITH  FIGURE  8 SHOWING  THE  SAME  RUN 
WITH  TRANSITION  DETERMINED  BY  THE  BOUNDARY  LAYER 
PROBE 


FIG.  8 SAME  RUN  AS  IN  FIG.  7.  A COMPARISON  OF  FIGS.  7 S 8 
SHOWS  THAT  THE  BOUNDARY  LAYER  PROBE  INDICATED 
TRANSITION  AT  THE  SAME  R*x  AS  THE  FRICTION  FACTOR 
MEASUREMENTS. 


LOCAL  APPARENT  FRICTION  FACTOR,  4fApp  LOCAL  APPARENT  FRICTION  FACTOR,  4f 


FIG.  9 ADIABATIC  PERFORMANCE  OF  THE  SYSTEM  FOR  SEVERAL 
FLOW  RATES  AND  RATES  OF  BOUNDARY  LAYER  SUCTION. 


LOCAL  APPARENT  FRICTION  FACTOR,  4fApp  LOCAL  APPARENT  FRICTION  FACTOR,  4f 


FIG.  II  THE  EFFECT  OF  HEATING  ON  TRANSITION.  A COMPARISON 
WITH  FIG.  10  INDICATES  A DECREASE  IN  TRANSITION  Rex 
DUE  TO  ACCUMULATION  OF  DIRT  PARTICLES  ON  WALLS. 


FIG.  12  THE  EFFECT  OF  HEATING  ON  TRANSITION.  SAME  AS 
FIG.  II,  EXCEPT  THAT  MEAN  BOUNDARY  LAYER  REYNOLDS 
NUMBER  IS  USED  IN  PLACE  OF  FREE  STREAM 
REYNOLDS  NUMBER. 


FIG.  13  ADIABATIC  RUN  USING  BOUNDARY  LAYER  PROBE.  RUN 
TAKEN  AFTER  SYSTEM  HAD  JUST  BEEN  THOROUGHLY 
CLEANED. 


FIG.  14  THE  EFFECT  OF  BOUNDARY  LAYER  SUCTION  RATE  ON 
THE  ADIABATIC  PERFORMANCE  OF  THE  SYSTEM.. 
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FREE  STREAM  REYNOLDS  NUMBER,  Re* 


FIS.  15  THE  EFFECT  OF  HEATING  ON  TRANSITION  FOR  A FLOW 
RATE  GIVING  R«0  * 130,000- 


FIG.  16  THE  EFFECT  OF  HEATING  ON  TRANSITION.  SAME  AS 
FIG.  15,  EXCEPT  THAT  MEAN  BOUNDARY  LAYER  REYNOLDS 
NUMBER  IS  USED  IN  PLACE  OF  FREE  STREAM 
REYNOLDS  NUMBER. 


FIG.  17  THE  EFFECT  OF  HEATING  ON  TRANSITION.  THE  EFFECT 
OF  BOUNDARY  LAYER  SUCTION  RATE  ON  THE  AOIABATIC 
PERFORMANCE  OF  THE  SYSTEM  FOR  THE  SAME  R«0 
IS  ALSO  SHOWN. 


FIG.  18  THE  EFFECT  OF  HEATING  ON  TRANSITION.  SAME  AS 
FIG.  17,  EXCEPT  THAT  THE  MEAN  BOUNDARY  LAYER 
REYNOLDS  NUMBER  IS  USED  IN  PLACE  OF  THE 
FREE  STREAM  REYNOLDS  NUMBER. 


TOTAL  rSTATIC  • 


FIG.  19  THE  EFFECT  OF  HEATING  ON  TRANSITION  FOR  A LOW 

FLOW  RATE  WITH  R«0  * 30,000.  THE  EFFECT  OF 

BOUNDARY  LAYER  SUCTION  RATE  ON  THE  ADIABATIC 

PERFORMANCE  OF  THE  SYSTEM  FOR  THIS  Rt0  IS 

ALSO  SHOWN. 


MEAN  BOUNDARY  LAYER  REYNOLDS  NUMBER,  RtXm 


FIG.  20  THE  EFFECT  OF  HEATING  ON  TRANSITION.  SAME  AS 
FIG.  19,  EXCEPT  THAT  THE  MEAN  BOUNDARY  LAYER 
REYNOLDS  NUMBER  IS  USED  IN  PLACE  OF  THE 
FREE  STREAM  REYNOLDS  NUMBER. 


FIG.  21  THE  EFFECT  OF  HEATING  ON  TRANSITION.  CLEANING  THE 
SYSTEM  BEFORE  THE  RUNS  WERE  TAKEN  HAS  CAUSED 
THE  VALUE  OF  R«x  AT  TRANSITION  TO  BE  TWICE  THE 
VALUE  SHOWN  IN  FIGS.  19  ft  20. 


FIG.  22  THE  EFFECT  OF  HEATING  ON  TRANSITION.  SAME  AS 
FIG.  21,  EXCEPT  THAT  THE  MEAN  BOUNDARY  LAYER 
REYNOLDS  NUMBER  IS  USED  IN  PLACE  OF  THE  FREE 
STREAM  REYNOLDS  NUMBER. 


TOTAL  ^STATIC  • 


FI6. 23  THE  EFFECT  OF  HEATING  ON  TRANSITION.  A VERY 
THOROUGH  CLEANING  HAS  INCREASEO  THE  ADIABATIC 
R*x  AT  TRANSITION  TO  850,000. 


FIG.  24  THE  EFFECT  OF  HEATING  ON  TRANSITION.  SAME  AS 
FIG.  23,  EXCEPT  THAT  THE  MEAN  BOUNDARY  LAYER 
REYNOLDS  NUMBER  IS  USED  IN  PLACE  OF  THE 
FREE  STREAM  REYNOLDS  NUMBER. 
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FREE  STREAM  REYNOLDS  NUMBER,  Rtx 


FIG.  25  THE  EFFECT  OF  HEATING  ON  TRANSITION.  INCREASING 

THE  TEMPERATURE  DIFFERENCE  BETWEEN  THE  WALL  AND 
FREE  STREAM  HAS  CAUSED  A CORRESPONDING  INCREASE 
IN  THE  UPWARD  DISPLACEMENT  OF  THE  CURVES. 


MEAN  BOUNDARY  LAYER  REYNOLDS  NUMBER,  R«*m 


FIG.  26  THE  EFFECT  OF  HEATING  ON  TRANSITION.  SAME  AS 
FIG.  25,  EXCEPT  THAT  THE  MEAN  BOUNDARY  LAYER 
REYNOLDS  NUMBER  IS  USED  IN  PLACE  OF  THE  FREE 
STREAM  REYNOLDS  NUMBER.  INCREASING  THE  TEMPERATURE 
DIFFERENCE  HAS  CAUSED  AN  INCREASE  IN  THE  MEAN 
BOUNDARY  LAYER  LENGTH  REYNOLDS  NUMBER  OF  TRANSITION. 


FIG.  27  THE  EFFECT  OF  HEATING  ON  TRANSITION.  THESE  CURVES, 
TAKEN  AT  A LOWER  FLOW  RATE,  EXHIBIT  THE  SAME 
CHARACTERISTICS  SHOWN  IN  FIG.  25. 


FIG.  28  THE  EFFECT  OF  HEATING  ON  TRANSITION.  SAME  AS 
FIG.  27,  EXCEPT  THAT  THE  MEAN  BOUNOARY  LAYER 
REYNOLDS  IS  USED  IN  PLACE  OF  THE  FREE  STREAM 
REYNOLOS  NUMBER. 


FIG.  29  ADIABATIC  CURVES  SHOWING  THAT  THE  ADDITION  OF 
DETERGENT  TO  THE  SYSTEM  DID  NOT  INFLUENCE 
TRANSITION. 
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FIG.  30  COMPARISON  OF  THE  RESULTS  OF  THEORETICAL  FRICTION 
FACTOR  COMPUTATIONS  FOR  THE  LAMINAR  ENTRY  OF  A 
TUBE  WITH  THE  DATA  OF  KUNE  ANO  SHAPIR0(6). 


FIG.  32  GENERALIZED  CURVE  OF  LOCAL  APPARENT  FRICTION  FACTOR  IN  A LAMINAR 
TUBE  ENTRY 


FLAT  PLATE : 

POHLHAUSEN  FLAT  PLATE 
VELOCITY  PROFILE 


FIG.  34  BOUNDARY  LAYER  GROWTH  IN  A LAMINAR  TUBE  ENTRY  AS  COMPARED  WITH 
THE  GROWTH  ON  A FLAT  PLATE 


